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ABSTRACT 

This study develops cost estimating relationships (CER) to 
estimate the cost of supply support for the maintenance of 
aircraft simulators. These CER's will be used in the 
Contractor Operation and Maintenance of Simulators (COMS) 
environment. The analysis techniques used to develop and 
validate these relationships include linear regression (simple 
and multiple), analysis of univariate distributions, 
simulation modeling, and linear optimization modeling. The 
regression analysis concludes that no useful CER is present in 
the cost elements represented by the sample of data gathered 
from actual simulator operations. There is, however, a useful 
CER present in another smaller set of data derived from the 
successful bids of eight contractors. These submissions were 
for maintenance with and without supply support. A simulation 
model was constructed to provide an independent cost estimate 
for use with the maintenance CER developed above. Analysis of 
univariate distributions was used to transform maintenance 
data from an operational simulator suite for use in the 
verification of the simulation model. Finally, an 
optimization model formulation is recommended for further 
investigation to determine the best mix of contractor-provided 
and Navy-provided spare parts to complete the optimization of 


Navy expenditures for supply support for aviation simulators. 





I. INTRODUCTION 


A. BACKGROUND 

In 1969 all aviation simulators were maintained and 
operated by personnel with a Tradevman (TD) rating. The Navy 
supply system was responsible for providing all repair parts 
and replacement assemblies for simulators. This system 
provided an adeguate number of operational ready (OR) 
simulator hours to meet the needs of the various aviation 
communities throughout the Navy. 

In the early 1980's, however, the Department of the Navy 
dramatically changed the way aviation simulators were 
Maintained and operated. The Defense Officer Personnel 
Management Act (DOPMA), force level constraints and the goal 
to increase the size of the Navy to 600 ships were the major 
factors that prompted the changes. In order to man 600 ships, 
many new “at sea" billets would be needed, and the Navy end 
strength was already at the maximum level permitted by law. 
Decisions by the Congress forced these additional billets to 
be generated internally. The additional "at sea" manpower was 
generated by replacing much of the military-manned permanent 
shore services with civilian commercial contracted services. 

The maintenance and operation of simulators became a prime 


target for conversion to commercial contract, and the concept 


of Contractor Operation and Maintenance of Simulators (COMS) 
was born. 

The conversion to COMS started in earnest in 1984 and was 
virtually complete by 1985 when the TD rating was 
disestablished. The original contracts were for maintenance 
and operation only, with the Navy continuing to be responsible 
for supply support. The COMS contracts were written so that 
the contractor was not charged down time when the Navy was not 
able to supply the reguired repair parts. 

The maintenance of a complete stock of repair parts is 
difficult for the Naval supply system. The large number of 
parts unique to simulators, coupled with their low utilization 
rates, places an undue burden on a system which is not 
designed to handle parts on a small scale. Тһе DOD 
acquisition regulations further hamper efforts to have 
sufficient spare parts to support the COMS contractor. The 
requirement for competition, or waiver for sole source 
acguisition, and the lead time reguired by the contracting 
process itself, freguently cause unusually high prices and 
unacceptable delays in the arrival of vital repair parts. The 
effect of these delays can cause an excessive number of 
unusable trainers and lead to decreased readiness for Naval 
aviation units. 

One solution to this dilemma may be a new generation COMS 
contract that makes the contractor responsible for some level 


of supply support. To aid in the conduct of a feasibility 


study to support the implementation of this new COMS contract, 
it is necessary to devise methods to answer the following 
guestions: 

1. How much does it cost to provide supply support? 

2. What is the cost of a COMS contract? 

3. Is there a best mix of contractor-supplied or 

Navy-supplied spare parts? 
B. SCOPE AND LIMITATIONS 
1. Scope 

This research focused on the development of two CERs 
that can be used to estimate the cost of spare parts and the 
cost of a COMS contract required to maintain a simulator for 
one year. 

This thesis summarizes the analysis and presents the 
results. Additional ideas are presented that outline the 
procedures needed to determine the best mix of Navy 
furnished/contractor furnished spare parts. 

2. Limitations 

The accuracy of the results of this analytical process 
is no better than that of the data used. Further limitations 
exist because of the small number of available simulators with 
useful data for this project. Other specific limitations will 
be discussed throughout this thesis as they occur. 

The intent of this thesis is to demonstrate analytical 
procedures that lead to the determination of a CER. The 


accuracy of the resulting CERs is not expected to be 


sufficient to allow its use in the contracting or budgeting 


process. 


II. COST ESTIMATION 


A cost estimate is a judgment or opinion regarding the 
cost of an object, commodity, or service. This judgment or 
opinion may be arrived at formally or informally by a 
variety of methods, all of which are based on the 
assumption that experience is a reliable guide to the 
future." [Ref.1:p.1] 

Cost estimation techniques are diverse and range from 
intuition to detailed application of labor and engineering 
cost standards. The five techniques most prevalent in The 
Department of Defense are industrial engineering standards; 
rates, factors, and catalog prices; cost estimating 
relationships; specific analogies; and expert opinion. There 
are many more technigues in the civilian sector. The 
determination of the technigue to be applied is dependent on 
the data available, what the estimate is to be used for, the 
time available to do the study, and other factors that affect 
the direction of the study. For the purposes of this thesis, 


cost estimating relationships and industrial engineering 


standards will be the only two methods explored. 


A. COST ESTIMATING RELATIONSHIPS 

Cost estimating relationships (CERs) are equations derived 
from statistical analysis that are used to estimate the cost 
of a product based on some characteristic of that product. 
This method is preferred because of its ease of application 


and accuracy. Its use requires the availability of high 


guality data and the existence of relationships among elements 
of that data base. A CER can be generated using regression 
analysis processes. The details of that process are discussed 


in the next chapter. 


B. INDUSTRIAL ENGINEERING STANDARDS 

On some occassions sufficient data are not available, or 
a useful CER can not be derived from the data, or a change in 
the process or materials makes historical data unusable. In 
these instances industrial engineering standards must be used. 

Estimating by industrial engineering procedures can be 
broadly defined as an examination of separate segments of 
work at a low level of detail and a synthesis of the many 
detailed estimates into a total. [Ref. 1:p.2] 

In this estimation scheme the process of production is 
broken down into its smallest elements. The cost of labor and 
material for each of these elements of production is then 
computed and totaled. The results are then multiplied by 
factors, determined by management, to estimate overhead and 


other indirect costs. The summation of these direct and 


indirect costs equals the total cost estimate for the product. 


C. PROBLEMS 

As is the case with all estimation processes, there are 
some inherent problems in each of these procedures. In the 
case of cost estimating relationships, there are two 
limitations. The first is the inherent uncertainty always 


associated with the use of statistics. The second is the 


uncertainty of the application of a particular CER to a 
certain situation. 

The industrial engineering standards technique has three 
drawbacks. First, it takes a large number of personnel with 
a detailed knowledge of the production process to compute the 
production element costs. Second, small mistakes in the 
production element costs can be magnified when multiplied by 
the predetermined factors in the computation of indirect 
costs. Third, the whole cost frequently turns out to be 
greater than the sum of all the production element costs. 

In considering the benefit to be derived from a cost 
estimation, one must always keep in mind that the estimations 
must be reasonable and structurally sound. The usefulness of 
the cost estimation may also be affected by the analyst 
himself. The lack of personal knowledge or familiarity with 
the systems being considered or a personal bias interjected 
into the analysis by the analyst can detract from the 


credibility of the estimate. 


III. METHODS OF ANALYSIS 


A. DETERMINING COST ESTIMATING RELATIONSHIPS 
1. Simple Regression Analysis 

Simple regression analysis is defined and explained in 
books that discuss data analysis. In Bachelder's An 
Introduction to Equipment Cost Estimating [Ref. 1], simple 
linear regression is defined aS a simple two-variable model 
which describes the linear relationship between the two 
variables by the equation 

у = а + bx + e (321) 

where y represents the dependent variable (cost in this case), 
x represents the independent variable (any of the other data 
elements), and e is a random error or fluctuation. The 
constants a and b represent the Y-axis intercept and the slope 
of the line, respectively. Eguation (3.1) represents the 
linear relationship between x and y, and it can be used to 
estimate y if an accurate value is available for x. The 
parameters a and b are estimated from data pairs (х;, 
Yi)... (X„, Yn) using well established analytical methods. The 
equations for these estimates A and B for a and b, 


respectively, are as follows: 


when = ір, =- Жу: 
П 
The STSC program Statgraphics was used in this thesis to 
perform the statistical analyses. This program uses the 
ordinary least sguares method to estimate the values of a and 
b, and will perform the calculations need to make statistical 


inferences about the parameters a and b and the line a + bx. 


a. Standard Error of the Regression (SE) 
The standard error of the regression is a value 
that represents the standard deviation of the estimated values 
of the error, e. The equation for SE in terms of the data 


pairs (X;, y;) is 


SE=/2 (y,- (A+Bx,))? 


where A and B are the least sguares estimate of a and b, 


respectively. 


b. T-value 
The t-value associated with the parameter b is 


defined as 


> دم 
SE [E (EX).‏ 


It is used to test the significance of the relationship 
between x and y. This is done by performing a statistical 
test of the hypothesis that the slope coefficient, b, is equal 
to 0. That is, it is decided that ЬЫ # O when the value of 
the t-statistic in equation 3.2 is larger in absolute value 


than the value of a number, t taken from standard 


а, п-2! 

student-t tables. The value of a is the level of significance 
for this test of hypothsis. 

с. R-squared Value 

The R-squared value is a measure of how well the 

regression equation predicts the value of y. It is frequently 

called the coefficient of determination, and equals one if all 

of the observed y, values are on the predicted line of the 

estimated model, A + Bx. 

2. Multiple Regression 

If the regression of one variable on another fails to 

produce an adequate CER, then a technique called multiple 

regression can be used. This method is similar to simple 

regression except it has more than one independent variable. 


In the resulting representative equation 


10 


y DRE CI eS, (OO) 
y is still the dependent variable, x and z are independent 
variables, e is the random error or fluctuation, a is the 
Y-axis intercept, and the coefficients b and c represent the 
change in y for each change of one unit of their respective 
variables (x and Z). 

It is important to remember that each added variable 
reduces the accuracy of the statistical inferences made about 
the paramters a, b, c and about the line itself. This becomes 
critical when working with small data sets as it can decrease 
the credibility of the analysis. 

The computed t-value and R-squared numbers are used 
for the same purposes as in simple regression. 

a. Collinearity 

Collinearity exists when two of the variables in 
a multiple regression model have, between themselves, a linear 
relationship. An indication that this condition exists is the 
presence of a high R-squared value while one or more of the 


variables has a low t-value. 


B. SIMULATION MODELING TECHNIQUES 

If statistical analysis fails, another form of analysis 
known as Simulation modeling may be used. Simulation modeling 
is the second technique to be utilized in this thesis. 

One dictionary definition of the verb simulate is as 
follows: "to assume or have the appearance or characteristics 
IR SN 


T 


This would appear to be a reasonably accurate description of 
the technique. The key to a successful simulation, however, 
is not just the "appearance or characteristics," but rather 
the operation or behavior of the simulation model. 

A useful working definition of simulation is the 
following: 

The process of designing a mathematical or logical model 
of a real system and then conducting computer-based 
experiments with the model to describe, explain, and 
predict the behavior of the real system. [Ref. 3:p. 5] 

Utilization of this second technique requires the 
construction of a model of the aircraft simulators that 
Simulates the maintenance process. This model utilizes 
parameters, unique to the performance of the system being 
investigated, to estimate the number of man-years required to 
Maintain the system. This information is then used to 
determine the cost of a COMS contract for a year. 

To validate this model, it is necessary to analyze 
historical data from a known system and use the analysis 
results to run the simulation. The results of the simulation 
are compared to reality. If the two are statistically 
equivalent, the model is validated and usable. 

1. Analysis of Maintenance Data 

The analysis of the maintenance data is fairly 
straightforward. First, it is necessary to obtain a data base 
that contains the time each failure occurred and the 
maintenance time it took to repair each discrepancy. The time 
of failure data is used to determine the inter-arrival times 
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(time between failures). The inter-arrival times and repair 
times are then analyzed to determine their probability 
distributions and the parameters of those distributions. 
These parameters are fed into the model which then simulates 
the maintenance process. 
a. Histogram 

The first step in analysis is to display a 
histogram of the data. A histogram is a graph that shows the 
number of occurrences of data points within a bin (range of 
values). The histogram will indicate the location (mean) and 
Variance of the probability distribution associated with the 
data. It may also give additional indications about this 
probability distribution. 

b. Probability Plot 

The probability plot is the primary means of 
determining the distribution reflected in the data set. This 
plot compares the quantile of a data set with the theoretical 
quantile of a selected distribution. If the two plots are 
reasonably close, the probability distribution that generated 
the sample is said to be the same, statistically, as the 


theoretical distribution. 


T3 


c. Kolmogorov - Smirnov Test 
The Kolmogorov - Smirnov test is used to determine 
the goodness-of-fit of a sample distribution to a theoretical 
distribution. It does this by comparing the absolute maximum 
distance between the sample cumulative distribution function 
and the theoretical cumulative distribution function. The 
computed significance level is used to determine if the sample 


distribution is the same as the theoretical distribution. 
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IV. COST ESTIMATION MODELS 


A. COST DATA 

The data used in this thesis were made available by Fleet 
Aviation Specialized Operational Training Group, Atlantic 
Fleet (FASO). The data represent the annual costs incurred in 
the acguisition and operation of forty-six separate aviation 
simulators for FY-87. The data elements provided include the 
cost of device acguisition (ACOCOST), number of hours the 
simulator was programmed for use (PGMHRS), COMS contract cost 
(COMCOST), utilities cost (UTILCOST), spare parts cost 
(SUPPCOST), total cost (TOTLCOST), and the cost per hour to 
operate the simulator (HRLYCOST). TOTLCOST is the sum of 
COMCOST, UTILCOST, SUPPCOST, and overhead costs for contract 
compliance inspection personnel. HRLYCOST is TOTLCOST divided 
by PGMHRS. 

Additional data on the cost of eight five-year COMS 
contracts bid with and without supply support were provided by 
Naval Training Systems Center (NTSC), Orlando, Florida. 

The first step in the analysis of this data was to see if 


a CER existed that would estimate the cost of supply support. 
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B. DEVELOPMENT OF CER USING FASO DATA 
1. Simple Regression 

Using the cost data supplied by FASO, a simple 
regression of SUPPCOST against all the other data elements 
demonstrated the absence of a useful CER with any of them. 
For the purposes of comparison, the t-statistic for 45 degrees 
of freedom and level of significance a=.05 is 2.414. A 
statistical summary of these regressions is provided in TABLE 
4.1. Note that the t-value of the slope in all cases but 
PGMHRS is larger than the t-statistic indicating that they are 
all significant, implying that the b is not egual to O. 
However, the R-squared values are all much lower than that 
required for a reasonable explanation of the dependent 
variable. 


TABLE 4.1 SIMPLE REGRESSION OF INDEPENDENT VARIABLES AGAINST 
THE DEPENDENT VARIABLE SUPPCOST 


INDEPENDENT PARAMETER T R 

VARIABLE VALUE SOUARED 
UTILCOST SLOPE 4.88562 35. 17% 
HRLYCOST SLOPE 7.68515 57.38 
TOTLCOST SLOPE 9.61059 67.735 
ACOCOST SLOPE 4.94474 35:72 
PGMHRS SLOPE . 42723 „41% 
COMCOST SLOPE 3.60971) 22.84% 


2. Multiple Regression 
With the failure of simple regression to produce a 


reasonably accurate CER, a backward stepwise multiple 
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regression technigue was tried. Since SUPPCOST is a component 
of both TOTLCOST and HRLCOST, it would have a linear 
relationship with both of those variables. Therefore, in 
order to insure that collinearity was not a problem, TOTLCOST 
and HRLCOST were not considered during the multiple 
regressions. 

TABLE 4.2 shows the results of the multiple regression of 
SUPPCOST against the remaining four variables. The first step 
yields little significance for any of the variables, and the 
R-sguared value is very low at 34.21. The next step 
eliminates the variable with the lowest t-value (COMCOST) and 
does the regression process again. The results of this step 
are not much better. The significance of each of the 
variables is less, and the R-sguared value does not change 
appreciably. These results indicate that these data do not 
exhibit the relationships reguired to produce a CER that can 
be used to estimate spare parts cost. 

The failure of the data to produce a usable CER is not 
surprising. The funding for spare parts comes from 
operational funds. For this reason, the cost of supplies is 
a measure of the availability of money as well as the cost of 


replacement repair parts. 
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TABLE 4.2. MULTIPLE REGRESSION OF INDEPENDENT VARIABLES 
AGAINST THE DEPENDENT VARIABLE SUPPCOST 


RUN 1 
Model fitting results for: SUPPCOST 

Independent 

variable coefficient std. error t-value sig.level 
CONSTANT 28435.122998 2.039639Е41 1.3941 0.1708 
PGMHRS -8.033911 6.666066 =L. 2052 0.2350 
ACOCOST 0.002245 0.001477 J 5196 0. T363 
COMCOST -0.049499 0.110375 -0.4485 0.6562 
UTILCOST 0.270332 0.205072 1 3182 0.1947 


R-SQ. (ADJ.) = 0.3421 SE= 38015.056957 MAE= 20951.79413 
Previously: 0.0000 0.000000 0.000000 
46 observations fitted. 


RUN 2 

Model fitting results for: SUPPCOST 
Independent 
variable coefficient std. error t-value sig.level 
CONSTANT 25588.820008 1.919855E4 1 9 0.1898 
PGMHRS 7.789751 6.580307 =T; 1838 0.2432 
ACQCOST 0.002071 0.001412 1.4669 0.1499 
UTILCOST 0.232652 0.185282 12557 0.2162 


R-SQ. (ADJ.) = 0.3546 SE= 37651.779438 MAE= 20686.445263 
Previously: 0.3421 38015.056957 20951.794135 
46 observations fitted. 
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C. DEVELOPMENT OF CER USING NTSC DATA 

The data provided by Naval Training Systems Center (NTSC) 
consists of the successful bid for each of eight different 
simulator maintenance contracts. These bids included the cost 
of five years of maintenance both with and without supply 
support making it most applicable to the questions at hand. 

The results of a simple regression of COMS cost without 
Supply support (COMWOSUP) as the independent variable and cost 
of supply support (SUPPLY) as the dependent variable were very 
good. The t-value of the slope was 9.88510, compared to a 
tabled value of 2.998 for 7 degrees of freedom with an 
alpha=.01, indicating the significance of the coefficient. 
The R-squared value 1s a respectable 94.21. This is probably 
not good enough for budgeting purposes, but it is accurate 
enough to be used for the comparison of contracting 


strategies. The resulting CER is: 
Cost of Supply Support= -37012.8 +.1565*COMS Cost +e (4.1) 


It should be noted that the cost for supply support and the 
COMS cost are for a full five year contract. 

The problem with this CER is that it is now necessary to 
come up with a method to determine an independent estimate for 


the cost of a COMS contract before it can be useful. 
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D. DEVELOPMENT OF CER FOR COMS CONTRACT COST 

The regression analysis procedures were repeated with 
COMCOST as the dependent variable. The results were similar 
to the results of the regression with SUPPCOST above. The t- 
values were good for some of the regressions, but the R- 
sguared values were not nearly high enough to indicate a 
guality CER. The multiple regression for COMCOST was not 
successful either. Stepwise regression also produced low R- 
sguared values, again indicating a less than desirable CER. 
Hence, no CER exists in the FASO data set for the cost of a 


COMS contract. 


E. SUMMARY OF REGRESSION RESULTS 

Regression has yielded a CER for the cost of supply 
support. That CER uses the cost of the COMS contract as an 
independent variable. Regression failed to produce a usable 
CER for the cost of a COMS contract, so another method must be 
utilized to provide that cost. In the next chapter, such an 


alternate procedure will be demonstrated. 
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V. SIMULATION MODEL 


The failure of the previous data sets to produce a CER 
provides an opportunity to explore the modeling of engineering 
standards to provide a method of cost estimation. 


An old proverb states that a journey of a thousands miles 
begins with but a single step. To paraphrase, a simulation 
analysis begins with a single step and must proceed 
methodically through several steps if useful results are to 
be derived and implemented. But before beginning any 
significant journey, it is important to know something about 
the route and noteworthy milestones along the way. 


Simulation analysis is a descriptive modeling technigue. As 
such, it does not provide the explicit problem formulation 
and solution steps which are an integral part of 
optimization models, such as linear programming. 
Consequently, one must specify in some detail a procedure 
for the development and use of simulation models to assure 
successful outcomes from their application. [Ref. 3:pp. 14- 
15] 


The process utilized in this thesis to insure accurate 
results consists of 

1) A definition of the problem and understanding of the 

process to be simulated so a good model formulation can 


be constructed; 


2) A search for data and analysis of that data to determine 
its characteristics; 


3) Development of a model to simulate the maintenance 
process; 


4) Verification and validation of the model's code, logic 
and results; 


5) Utilizing simulation results to determine the cost of a 
Simulator maintenance contract. 
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A. PROBLEM DEFINITION AND FORMULATION 

Problem formulation is, in a sense, the most important 
step in a simulation analysis. Appropriate solutions to 
inappropriately formulated problems cannot be achieved. 
But before a problem can be formulated, it must be 
identified or found. Problem finding is, in reality, 
choosing from among several problems which are competing 
for the same resources. Criteria for selection include 
technical economic, and political feasibility, and the 
perceived urgency for a solution. Problem selection can 
have a significant impact on the ultimate success of the 
analysis and the implementation of the results. [Ref. 3:pp 
16-17] 

In order to formulate this simulation model, it is first 
necessary to have an understanding of the process to be 
Simulated and define the problem to be solved. 

1. The Simulator Maintenance Process 

There are two types of simulator maintenance, 
scheduled maintenance and nonscheduled maintenance. Scheduled 
maintenance is done as required by the preventive maintenance 
schedule, is not variable, and will therefore not be part of 
the model. 

Nonscheduled maintenance occurs as the result of a 
component failure or system breakdown. These component 
failures or system breakdowns can be discovered and reported 
in one of three scenarios. First, maintenance personnel can 
find a discrepancy during scheduled maintenance and report it. 
Second, a pilot can find a discrepancy during the conduct of 
a training mission and not report it until the end of the 
trainer period. Finally, a pilot can find a discrepancy 
during the conduct of a training mission that degrades the 


performance of the simulator below that reguired for the 
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completion of the mission, and the pilot reports that 
discrepancy immediately. All three of these scenarios are 
called breakdown events and will be assumed to occur randomly. 

After a discrepancy has been reported, it is assigned 
to a technician for repair. This repair takes a variable 
length of time to complete. The completion of the repair 
makes the simulator available for use and is called a repair 
event. The simulator is scheduled for use and the process 


starts over again. 


2. Problem Statement 
To answer the questions posed in this thesis, it is 
necessary to determine how many man hours are required to 
Maintain an aviation simulator. This will be converted to the 
number of man years required for maintenance, and used to 


determine a cost estimate. 


B. SIMULATION MODEL DEVELOPMENT 

The problem of determining the number of man hours 
required to maintain an aviation simulator can be answered by 
developing a Fortran language computer model to simulate the 
maintenance of an aviation simulator suite. This model should 
be easily modifiable to fit simulator suites consisting of any 
number of elements. An element is a segment of a simulator 
suite designed for a specific training mission (i.e., 


operational flight trainer (OFT), weapons system trainer 


(WST), etc). 
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Such a model was developed for this thesis. It utilizes 
a main program to initialize the parameters and perform the 
simulation by calling a variety of subroutines. These 
subroutines are the initialization subroutine, the time 
advance subroutine, the breakdown event subroutine, the repair 
event subroutine, the report generator subroutine, and random 
variable generators for the exponential, uniform, and normal 
distributions. The failure and repair subroutines are 
constructed from a block of code with one block for each 
Simulator element. This enables the model to be expanded or 
contracted to meet the specifications of any simulator suite. 
At the completion of the simulation, the report generator 
prints out the parameters used for the failure and repair 
distributions, the number of maintenance personnel, maximum 
queue length, element operational ready rate, number of 
maintenance actions that took longer than three hours to 
repair, Simulation run time, and total number of maintenance 
actions that occurred during the simulation. A complete 
listing of the program is included in Appendix A. 

1. How to Use the Model 


a. The first step is to load the program into a computer 
with a Fortran compiler. 


b. Modify the parameters for the failure and repair 
distributions to reflect those of the simulator suite 
being simulated. An analysis of these distributions for 
the example considered in this thesis is in the next 
section. 
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Modify the number of elements in the program to reflect 
the physical composition of the simulator suite of 
interest. For the purposes of this thesis, the eight 
elements, Operational Flight Trainer (OFT) 1 and 2, 
Weapons System Trainer (WST) 1 and 2, Positional Type 
Trainer Tactics (PTT) 1 and 2, and Positional Type 
Trainer Cockpit (PTC) 1 and 2 have been used. 


Estimate the number of maintenance personnel reguired to 
maintain the simulator suite, and enter that number for 
the value of NMTPER. Then, determine the length of time 
the simulation is to cover and enter that value for 
FPEL(17) 


Complle and run the model. 


Note the results of interest, in this case operational 
ready rate and maximum repair queue length, and modify 
the number of maintenance personnel up or down until the 
desired results are obtained. 


Determine the number of maintenance personnel reguired 
to perform the contract. This is done by multiplying 
the number of maintenance personnel needed to meet 
mission reguirements by 3 to cover a 3 shift 24 hour 
day, and adding the appropriate number of overhead 
personnel. Overhead personnel include, but are not 
limited to, site supervisor, supply clerk, maintenance 
administration clerk, personnel to perform the reguired 
preventive maintenance (PM), and safety observers as 
reguired by Navy safety regulations. 


Multiply the number of man-years reguired to maintain 
the simulator by the NTSC-provided cost for acguiring 
one man-year on contract, $45,000, to get a cost 
estimate for the COMS contract. 


DATA SEARCH AND ANALYSIS 


Data used in determining the parameters to be utilized by 


this model were obtained from an actual P-3 simulator site. 


Failure rate data were obtained for a 227 day period of time, 


and time to repair data were obtained for a different 90 day 


period of time. The different periods of time for the two 


data sets were driven by the way the records were kept by the 
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maintenance contractor. These data were analyzed to determine 
if the failure rate and repair times demonstrated the behavior 
of a known probability distribution. 
1. Analysis of Data 
a. Interarrival Time 

The periodic occurrence of an event such as the 
breakdown of an aviation simulator is best described in terms 
of the time between breakdown events. This period of time is 
called the interarrival time. The data used to determine the 
interarrival times were taken from the maintenance action 
forms (MAF) that were filled out by the person discovering the 
malfunction and submitted to maintenance personnel. Once the 
failure times were entered into a data bank, they were 
processed by a simple Fortran program to convert them into 
interarrival times. 

The analysis of the resultant interarrival times 
was done utilizing the IBM program GRAFSTAT. The "Fitting 
Probability Distributions" option of the “Analysis of 
Univariate Distributions" section of the Grafstat menu was 
used. This option utilized the fitted histogram/density and 
probability plots to analyze the data and determine its 
parameters. It was initially anticipated that the 
interarrival times would be exponentially distributed. 

The data used for this thesis turned out to be 
very ill behaved. The data did, in fact, fail to conform to 


any known distribution, and it defied initial attempts at 
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analysis. It was therefore necessary to attempt data 
transformations in order to get any usable information. 
(1) Transformation. Transformation is the 


systematic modification of data to make it conform to a 


distribution. The most utilized of these modification 
procedures are the logarithmic and exponential 
transformations. The process is a simple one; the data is 


raised to a certain power, or its logarithm is determined, and 
the resulting data set is analyzed. The exponential 
transformation proved to be useful and is represented by 
equation 5.1, 

y? = x. (5.1) 
Y is the data, a is the power to which the data is raised, 
and X is a random variable that can be fitted to some 
distribution. 

Unfortunately, even after transformation, much 
of the data continued to defy analysis. Probability plots of 
the raw data showed that much of it was bimodal or lacked good 
discipline. Neither of these problems came as a total 
surprise because of the inherent characteristics of the 
failure process. 

(2) Problems. A bimodal situation is created when 
data points on a data set come from two or more distributions. 
That is to say, in this case, that the failure rate of the 
main computer is different than that of the motion systen, 


which is different from that of the cockpit indications, etc. 
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When all the types of failures are thrown into one set of 
numbers, it can create many problems for the analyst. 

The data problems appear to be the result of 
the process by which the data were recorded. Many times a 
discrepancy will be discovered during a training evolution, 
but, since it does not degrade training, it is not recorded 
until the end of the exercise. Sometimes several 
discrepancies will be identified at essentially the same time 
(the completion of the trainer time period). There is alsoa 
high probability that several malfunctions were discovered by 
maintenance personnel and corrected without any paper work to 
record the malfunction or the repair. All of these types of 
data excursions cause the data to be ill behaved and increase 
the difficulty of analysis. 

(3) Solution. To counter these problems, several 
algorithms for the performance of the required analysis were 
devised and tested. The one that produced the best results is 
described as follows. 

The first step of this algorithm was to 
uniformly distribute all blocks of discrepancies obviously 
written at the end of a training evolution. This 
redistribution was made over the last half of a three hour 
training period. The resulting data were then analyzed as 
described above in two different ways. First, each element 


(OFT1, OFT2, etc.) was individually analyzed. Second, the 
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elements were combined by type (OFT, WST, etc.) and analyzed 
again. 

The results of this analysis are shown in 
TABLE 5-1. For each element and type, the second column shows 
exponent that was used to transform the interarrival times. 
The third column shows the value of the parameter of an 
exponential distribution which was fit to the transformed 
interarrival times for that element and type. The fourth 
column shows the significance level, using a Kolmogorov- 
Smirnov statistic, of the test of the hypothesis that the 
empirical distribution of the transformed interarrival times 
fit a theoretical exponential distribution with the parameter 


shown in the column three. 


TABLE 5.1 PARAMETRIC ANALYSIS OF THE INTERARRIVAL TIMES FOR 
SIMULATOR ELEMENTS 


Коло» 
ELEMENT TRANSFORMATION EXPONENTIAL SIGNIFICANCE 

EXPONENT PARAMETER LEVEL 
DET .84 21:59 9198 
52 40 4.9 3137 
OFT 88 18:92 < 7991 
WST1 . 64 1959 77 . 8061 
WST2 46 5767 . 9039 
WST 70 22.64 9057 
PTT1 54 8.86 . 5918 
PTT2 . 2 14745 6599 
PTT 7 28.05 9619 
PTICI 69 33465 6591 
PTC2 236 3.06 . 1296 
РТС 223 2199 „9256 
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The most significant results were then used in 
accordance with the following rules: 

1) If analysis resulted in good, significant 
parameters for both elements of a type of 
trainer, those individual parameters were used in 
the simulation. 

2) If analysis resulted in very significant results 
coming from only one element of a type, then that 
parameter was used for both elements in the 
simulation. 

3) If the most significant results came from the 
combined elements of type data, the parameter for 
the combined data was used, as is, for both 
elements in the simulation. 

In the analysis of the failure data, rule 2 
was applied to OFT1 and OFT2 (OFT1 data being used). Rule 3 
was applied to WST1, WST2, РТТ1,‚ PTT2, PTC1, and PTC2. The 
parameters used for interarrival times in the simulation are 
presented in TABLE 5.2. The effects of these choices of the 
parameters on the overall model will be shown in Section D 
below. 


TABLE 5.2 INTERARRIVAL TIME PARAMETERS USED IN THE 
SIMULATION MODEL 


ELEMENT TRANSFORMATION EXPONENTIAL 
EXPONENT PARAMETER 
ORTI . 84 9 15.9 
ОЕТ2 . 84 519 
WSTl ‚70 22.64 
WST2 210 22.64 
PEIL 77 28.05 
PEZ 77 28.05 
PTC 1 . 3 21.99 
PTC2 13 21799 
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b. Maintenance Repair Time 

The length of time it takes a maintenance 
technician to repair an aviation simulator system discrepancy 
is called the maintenance repair time. The maintenance repair 
times are taken directly from the portion of the MAF that is 
filled out by the technician after he has finished the repair 
process. 

The analysis of the maintenance repair time data 
was conducted identically to that of the interarrival time 
data, and the results were similar. The data turned out to be 
equally as ill behaved and the transformation process was 
again required to obtain usable results. 

The decision process was the same as well, and 
TABLE 5.3 illustrates the transformation exponents and mean 


service time parameters used in the model. 


TABLE 5.3 MEAN SERVICE TIME PARAMETERS USED IN THE 
SIMULATION MODEL 


ELEMENT TRANSFORMATION MEAN SERVICE STANDARD 

EXPONENT TIME DEVIATION 
OF'T1 33 1.14 2222 
OF TZ NATURAL LOG 2.08 . 36 
WST1 32 1.46 .43 
WST2 232 1.46 . 43 
ETTI "50 2423 . 74 
PTT2 250 2529 . 4 
P'TC I .14 1.05 .07 
PTC2 .14 1705 . 7 
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D. VERIFICATION AND VALIDATION OF THE MODEL 
Model verification and validation actually is concerned with 
three models: a conceptual model, a logical model, and a 
computer model. ...Verification focuses on internal 
consistency of a model, whole validation is concerned with 
the correspondence between the model and reality. (Ref. 3:p. 
2M 

Verification of the model was a continuing process 
throughout its development. As they were identified, code and 
logic errors were rectified until the model ran 
satisfactorily. 

To validate the simulation model, the parameters derived 
by the analysis of available maintenance data were inserted 
into the simulation model. The model was set to run for 6000 
hours which represents the number of hours of operation 
provided for in one year of operation (364 days - 104 weekend 
days - 10 govt holidays * 24 hours/day = 6000 hours). The 
program was run for two maintenance personnel (minimum number 
of personnel required by safety regulations for maintenance of 
energized equipment with high voltage components) and 
increased by one person until the maximum queue size was 0. 
The detailed results are shown in Appendix B. 

It appears that four maintenance personnel are required to 
keep the simulators up without any delay time due to non- 
availability of maintenance personnel. The preventive 
maintenance (PM) manuals indicate that it takes 3.15 people to 
do the required PM. 

All the information required to determine the number 


of maintenance personnel to maintain this simulator suit is 
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now known. That number is 18.15. TABLE 5.4 illustrates the 
determination of this number. 

The simulator in guestion was being maintained by 18 
people at the time the data was gathered. The model therefore 
appears to be validated as a good simulation of aviation 
Simulator maintenance. 


TABLE 5.4 NUMBER OF PERSONNEL REQUIRED TO MAINTAIN THE P-3 
SIMULATOR MODELED 


BILLET NUMBER REQUIRED 
SIGHT SUPERVISOR 1 
SUPPLY CLERK 1 
MAINTENANCE ADMINISTRATION CLERK l 
PM PERSONNEL 3s 15 
MAINTENANCE PERSONNEL 12 
TOTAL 18.415 


E. UTILIZATION OF MODEL RESULTS 

The solution to the problem of determining the cost of 
maintenance and supply support for aviation simulators is at 
hand. Simple computations are all that is reguired to 
estimate those costs. 

The cost of a contract reguiring 18.15 men to maintain 
simulators is computed by multiplying 18.15 by $45,000 (the 
approved estimate for the cost of one man year of an average 
simulator technician). That cost is $816,750 for one year of 
maintenance service. This is within 5$ of the cost of one 


year of service as delivered by the present contractor. 
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Now the developed COMS cost can be used in the developed 
CER (eguation 4.1) to determine the cost for contractor supply 
support for this simulator suite. This cost is estimated to 
be $127,089 per year. This also compares favorably with the 
cost of the contract now in effect. 

Therefore, one method of determining the cost of supply 


support provided by a COMS contractor has been determined. 
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VI. AREAS FOR FURTHER RESEARCH 


Presently the contractor supplied/Navy supplied mix for 
spare parts is determined by what is in the simulator supply 
room at the time the contractor takes over maintenance of the 
simulator. An inventory is conducted, and the contractor is 
responsible for maintenance of that inventory until the 
completion of the contract period. While this is a step in 
the right direction, there are ways of determining a more cost 
efficient spare parts mix that considers all non trainer 
peculiar parts. 

The best way to determine the most cost efficient parts 
mix is to construct a linear optimization program that 
Minimizes the cost of parts support. A proposed basic 


formulation for that optimization program is: 


INDEX 
1 = PART NUMBER 
PARAMETERS 
CNA = COST FOR NAVY ACQUISITION OF PART 
CNS = COST FOR NAVY TO STORE PART 
CNH = COST FOR NAVY TO HANDLE PART 
CCA = COST FOR CONTRACTOR ACQUISITION OF PART 
CCH = COST FOR CONTRACTOR TO HANDLE PART 


NR 


NUMBER OF PARTS REQUIRED FOR LIFE CYCLE OF 
SIMULATOR 
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VARIABLES 


X; = NUMBER OF PARTS , NAVY TO BUY 

Y; = NUMBER OF PARTS ; CONTRACTOR TO BUY 
FORMULATION 

MINIMIZE 


Z; CNA,*X, + CNS,*X, + CNH,*X, + CCA,*Y, + CCH,*Y, 


SUBJECT TO 
Хх + Y, = NR, 

This basic formulation is simplistic by design. It has no 
constraints, other than the parts constraint, to allow for 
tailoring to any specific situation. For example, if there is 
a limited amount of operating funds (O&MN) the following 
constraint can be added: 

2; CNA¡*X. < FUNDS AVAILABLE 

Likewise, if there is a paucity of funds available for 
contractor supply support, the following constraint can be 
added: 

2 CCA TY i CCH;*Y; < FUNDS AVAILABLE 

Values for X; can be assigned for those parts that are 
trainer peculiar and have to be purchased by the Navy at the 
time of acquisition of the simulator. 

The ability to tailor this formulation to any simulator 
suite is limited only by reality. If a constraint can be 


identified and linearly represented, it can be included in the 
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formulation. The only real limitations are the availability 
of data and the format of that data. 

Data were not available to support this model at the time 
the research was done for this thesis. The hard copy data 
that was acquired for a small subsystem did not have complete 
cost data for over-the-counter (contractor) purchase of most 
of the parts. A recommended format for the data to be 


optimized by this formulation is: 


PART COST 

INDEX ACQUISITION HANDLING STORE 
NUMBER NSN NAVY CNTR NAVY CNTR NAVY 
1 1234567890 145.00 128.00 1250 12750 13.40 


Once the data have been put into the right format, the 
formulation and its associated data can be an input into any 
matrix generator/linear solver. The resultant solution of the 
linear set of equations will provide the optimal mix of 
contractor/Navy acquisition of spare parts 

Unfortunately, the absence of usable data prevents the 
verification of this optimization model. However, the 
formulation's simplicity, coupled with its sound theoretical 
foundation, would lend an air of credibility to the results 


obtained. 
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VII. SUMMARY AND RECOMMENDATIONS 


A. SUMMARY 

The results of the analysis in this thesis show that the 
aircraft simulator cost elements studied do not correlate well 
enough to allow for successful linear regression analysis. 
This is most likely caused by pressure applied to the level of 
spending on repair parts from sources other than the 
maintenance process. Examples are reguirements to expend all 
funding during the quarter received and the necessity to 
transfer funds between commands to meet priority requirements 
in other areas. 

However, it has been shown that there is a predictable 
relationship between the cost of a maintenance contract and 
the cost for providing supply support for that contract. This 
relationship is derived from a relatively small set of data 
drawn from eight successful bidders for simulator maintenance 
contracts. The accuracy of this model is excellent, with a 
high R-squared value (94.21). Utilization of this model will 
provide reasonably accurate government estimates for 
comparison with the bids from commercial contractors. 

A maintenance process simulation model was created to 
provide an independent estimate for the cost of maintaining a 
Simulator for one year. Its accuracy was validated by 


analyzing the performance of an actual simulator and utilizing 
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those performance parameters in the simulation model. The 
results were impressive. In addition to its accuracy, the 
model, because of its modular construction, is easily 
adaptable to any specific simulator suite. 

This thesis has also provided a formulation for the 
determination of the optimal mix of contractor-supplied versus 
Navy-supplied repair parts. Even though it has not been 
verified or validated, this optimization formulation has great 
potential for savings if the tìme is taken to put the cost 


data in the format reguired. 


B. RECOMMENDATIONS 

The results of this study suggest that regression analysis 
and computer modeling can be of significant value in the 
estimation of supply support costs and maintenance costs for 
aircraft simulators. Continued utilization of these analysis 
techniques to update the CERs will further refine the accuracy 
of the results obtained. 

The potential cost savings that can be gained by 
utilization of the optimization formulation developed herein 
cannot be overemphasized. The fervent pursuit of a usable 
data base for all simulators in the inventory is highly 


recommended and will pay dividends in the future. 
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APPENDIX A 


& He TT e OAC e e e e TTI RIK IH IK RI ik K k جل جلا جلا‎ K K k k K k K kok k k جف جف جف جف جل جل جف جل جفة جف جف جل جف جلا‎ e de k kok k e e جف جف جف جف جف جف‎ De e de de e de de de de dr de de d de dee جف جف جف‎ Sie de de de die de de de K 
THSMOD 


A simulation model written to simulate the maintenance requirements for an aviation simulator suit 
and print out the results. 


WRITTEN BY 


CDR EDWARD M. BIGGERS 


TO SATISFY THE REQUIREMENTS FOR 


DEFINITION OF VARIABLES 


A MASTER OF SCIENCE DEGREE IN 
OPERATIONS RESEARCH 


SEPTEMBER 15, 1989 


BRKDN Break down subroutine 

CHKIN Time break down occurred 

CHKINO Time break down added to the gueue 
CLOCK Time keeper for simulation 

ELAWMT Simulator element waiting for maintenance 
EXPON Exponential random variate generator 
F Number of bread downs lasting longer that 3 hours 
FEL Future events list 

IAT Interarrival time 

IMEVT Imminent event 

INITLZ Initialize subroutine 

J,K,L Do-loop index 

MITOF 1 Mean interarrival time OFT1 

MITOF2 Mean interarrival time OFT2 

MITWS1 Mean interarrival time WST1 

MITWS2 Mean interarrival time WST2 

MITPT1 Mean interarrival time PTT1 

MITPT2 Mean interarrival time PTT2 

MITPC1 Mean interarrival time PTC1 

MITPC2 Mean interarrival time PTC2 

MQ Maximum length of queue 

MSTOF 1 Mean service time OFT! 

MSTOF2 Mean service time OFT2 

MSTWS1 Mean service time WST1 

MSTWS2 Mean service time WST2 

MSTPT1 Mean service time PTT! 

MSTPT2 Mean service time PTT2 


Wr rr k ŵr ŵr rrr rr K K k k ŵr k k K K K K K K K K K k Ve Ve de De de De De De De ŷr جل‎ Ve جل‎ Ve K K K K K K K K k ŵr ŷr rr ŷr K k r K K K K ŵr ŵr ŵr ŵr cr cr cr r K ŵr rir ŵr rr ŵr r r r r oir o r r r ŵr ŵr K K ŵr We ŵr ŵr Ŵr 
Wc جاو‎ dc e k kk k k k kk k kk kt Ve r rr ŵr KK k k Ve ŵr rr r r r r r K K K K Ve K k K k K K K K K K K ŵr e Ve Ve De cr cr cr cr r cr ŵr ŵr ŵr ŵr ŵr rr rr rr rr k k ŵr rr rr K K k k ŵr جل جل‎ Ne cr cr ŵr ŵr ŵr rr rr rr rc ŵr 


MSTPC1 Mean service time PTC1 

MSTPC2 Mean service time PTC2 

NAWTMT Number of simulator elements waiting maintenance 
NMTBUS Number of maintenance personnel busy 
NMTPER Number of maintenance personnel 
NORML Normal random variate generator 
NUMEVS Number of events 

OPRDY 1 Operational ready rate OFT1 

OPRDY2 Operational ready rate OFT2 

OPRDY3 Operational ready rate WST1 

OPRDY4 Operational ready rate WST2 
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OPRDY5 E Operational ready rate PTT1 


OPRDY6 = Operational ready rate PTT2 

OPRDY7 : Operational ready rate PTC1 

OPRDY8 E Operational ready rate PTC2 

REPAIR Е Repair subroutine 

RPTGEN Е Report generator subroutine 

SGMOF 1 - Standard deviation OFT1 service time 
SGMOF 2 - Standard deviation OFT2 service time 
SGMWS 1 - Standard deviation WST1 service time 
SGMWS2 - Standard deviation WST2 service time 
SGMPT 1 - Standard deviation PTT1 service time 
SGMPT2 - Standard deviation PTT2 service time 
SGMPC1 - Standard deviation PTC1 service time 
SGMPC2 - Standard deviation PTC2 service time 
SVT : Service time 

TIMADV : Time advance subroutine 

TNREP = Total number of repair events 

TREPTM : Total repair time 

TTMINQ š Total time in queue 


To utilize this model, first determine the parameters for each of the elements of the 
simulator. These parameters include interarrival time, mean service time for repairs, and the standard 
deviation for those repair times. Insert these values into the main program. Modify the break- 
down/repair subroutines to emulate the number of elements in the simulator suit being modeled. 
Determine the number of events and enter that number into the main program. 


Estimate the number of maintenance personnel required to maintain the simulator at the desired 
operational ready rate and enter that value in the main program. 


Run the simulation program and note the results. 


Modify the number of maintenance personnel until the operational ready rate is within an 


acceptable range, until the maximum queue size is 0, or until any other required condition is met. 
hr k k k k k k K k k k k k k k k k k k K k kt k k k k k k k k k k k k k k kt k k k k k k k k fr k k k k k fk k K k k ŷr Yr fk k k k k k k k k ŵr Yr ŵr ŵr k k k k k k k k k k k k k k k k k xr rr Wr vir її d k k x 
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PROGRAM THSMOD 
bc Se جل جلا جل جف جف جلا جلا جلا جلا جلا جف جف‎ e e جل جل‎ e e علا جلا علا علا جل جلا جل‎ 18+ OAC OAC OAC OAC OC ok OC OC O OC O ИКА He e e He Ae de de OAC de de He He e de He de e He de de e e e de de de e de e He Ae de He de K NA N NK AAA AAA AAI A 
MAIN PROGRAM 
1) INITIALIZE MODEL 
2) CALLS TIME-ADVANCE AND EVENT ROUTINES 
3) CALLS REPORT GENERATOR TO TERMINATE SIMULATION NORMALLY 


CREARA MARICA RIC ŵr de de ŵr ir rr de de irr rr Âr chr irr chr irr xr de de de de de rr cr rr rr chr rr cr rr rrr rr rr cr chr rr ŵr rr Âr dr rr âr rr rr rr rr kc rr ŵr ŵr rr ŵr rr ŵr x 


INTEGER F,ELAWMT(100),TNREP 
REAL MITOF1,MITOF2,MITWS1,MITWS2,MITPT1,MITPT2,MITPC1,MITPC2, 
1 MSTOF1,MSTOF2,MSTWS1,MSTWS2,MSTPT1,MSTPT2,MSTPC1,MSTPC2, 
2 FEL(17), CHKIN(8), CHKINQ( 100), NORML, TREPTM(8), TTMINQ(8) 
COMMON /SIM/ MITOF1,MITOF2,MITWS1,MITWS2,MITPT1,MITPT2,MITPCL, 
MITPC2,MSTOF1,MSTOF2,MSTWSL ,MSTWS2,MSTPT1,MSTPT2,MSTPC1,MSTPC2, 
SGMOF 1, SGMOF2, SGMWS1, SGMWS2, SGMPT1, SGMPT2, SGMPC1,SGMPC2,NAWTMT, 
NMTBUS, CHKIN, CHKINQ, TREPTM,MQ, F, TNREP, NMTPER, ELAWMT, TTMINQ 
COMMON /TIMEKP/ CLOCK, IMEVT,NUMEVS, FEL 
NUMEVS = 17 
Íc Ve dc de Ae de He He He de He He de de He He de de de e de de He he ke he he he h he ke k he he k k ke k he k k he k k k k k k k k k k k k k k k k k k Nr de de Be de de de de de de de de de de de Me Me Ne de Ve Ve LLL LLL LLL LLL LLL LL 
ASSIGN VALUES TO INPUT PARAMETERS 
(THESE VALUES COULD BE STORED IN A FILE AND READ INTO THE PROGRAM) 


ir cir ŵr ŵr ŵr Ae ŵr âr ŵr de r âr ŵr ŵr ŵr ŵr ŵr جل‎ ŵr جل‎ He ŵr ŵr de âc ŵr Âr ŵr ŵr de âc ŵr de ŵr K de de de Ae de de de ŵr He ŵr ŵr de De ŵr de kc ŵr ŵr ŵr de de de de De De De ŵr de de de De De de de De de de de de ŵr ŵr De de De De de he de c e eed de de 


UI N = 


ceo ee fef dede e e ee ee e ce e e ec ecce e e e ec ef ee e ee ecc K 0 k 0 RK * * * 1 * * * * * * * * * w * * * * * * * * * * * * * * * * * * * * * * * * * * 


MITOF1 = 25.799 
MITOF2 = 25.799 
MITWS1 = 22.64 
MITWS2 = 22.64 
MITPT1 = 28.053 
MITPT2 = 28.053 
MITPCL = 21.99 
MITPC2 = 21.99 
MSTOF1 = 1.14 
MSTOF2 = 2.0841 
MSTWS1 = 1.4604 
MSTWS2 = 1.4604 
MSTPT1 = 2.2272 
MSTPT2 = 2.2272 
MSTPCI = 1.051 
MSTPC2 = 1.051 
SGMOF1 = 0.2183 
SGMOF2 = 0.35516 
SGMWSt = 0.43201 
SGMWS2 = 0.43201 
SGMPT1 = 0.74425 
SGMPT2 = 0.74425 
ЅСМРСІ = 0.073798 
SGMPC2 = 0.073798 
NMTPER = 2 


CALL INITIALIZATION ROUTINE 


XRK ŵr ŵr rr ir ŵr rr ŵr ir ŵr جلا‎ ŵr ŵr ŵr ŵr ŵr ŵr ŵr rr rr ŵr rr ir dr 19: جل جل جل جل جل جل جل ل 117 117 جل جل جل جل جل جل 37 17 جل‎ ŵr rr âr ŵr âr rr r âr âr ŵr ŵr rr rr Âr de de de âr rr ŵr rr rr Wr ŵr ŵr âr ŵr rr ŵr ŵr ŵr ŵr ŵr rr ЖКК 


CALL INITLZ 
Tik kk k k k جف :ف‎ er sk جل جل ل جل جف‎ de + جل جف جا‎ K k ССС de de de de de de de de KOK OK OK ke ke ke 
CALL TIME-ADVANCE ROUTINE TO DETERMINE IMMINENT EVENT AND 
ADVANCE CLOCK TO THE IMMINENT EVENT TIME. 


re ŵr ir rr r rc Wr dcr rr rc dre ic ŵr ŵr ŵr rr ŵr ŵr ŵr âr ŵr ŵr rr rr rr de ŵr rr ŵr rr cr rr ccce eee e De r'r r'r r'r rr rr rrr ŵr rr rr rr rr rrr rr rr de de kc ŵr kc âr dee x 


30 CALL TIMADV 


XXX HX Á X КККК KK KK KKL 7ل‎ e ÉH ÉÉ K جل ف 117 جل 117 جل جل جل جل جل جل جل‎ 17 e e e e جل‎ x K 0 f K k É K k K k É * W * W É * * * * * * * * * * À À He À Àe Œœ * * * * 


THE VARIABLE 'IMEVT' INDICATES THE IMMINENT EVENT 


IMEVT = 1 IS AN OFT1 BREAK-DOWN 
IMEVT = 2 IS AN OFT2 BREAK-DOWN 
ІМЕУТ = 3 15 А WST1 BREAK-DOWN 
IMEVT = 4 ISA WST2 BREAK-DOWN 
IMEVT = 5 ISA PTT1 BREAK-DOWN 
IMEVT = 6 ISA PTT2 BREAK-DOWN 
IMEVT = 7 ISA PTCI BREAK-DOWN 
IMEVT = 8 ISA PTC2 BREAK -DOWN 
IMEVT = 9 IS AN OFT1 REPAIR 

IMEVT = 10 IS AN OFT2 REPAIR 

IMEVT = 11 ISA WST1 REPAIR 
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ІМЕУТ = 12 15 А WST2 REPAIR 
ІМЕМТ = 13 15 А РТТЇ REPAIR 
IMEVT = 14 ISA PTT2 REPAIR 
IMEVT = 15 ISA PTC1 REPAIR 
IMEVT = 16 ISA PTC2 REPAIR 


IMEVT 17 IS THE END OF THE SIMULATION 


r r k k k k'ik k k k k k k k k k k k kk k de de de K k k k k k K K kO k K de de k K k k O k k K KO KO R K de de K k K k de k k K k k k k k k k k k k k k k k k k k k k k k R k k k O RO RO KO KO ORO KO kO KO KO RO K K 


GO TO (40,40,40,40,40,40,40,40,50,50,50,50,50,50,50,50,60), IMEVT 


ŵr ŵr ŵr ŵr ŵr k'ik k k ik k k k k k k k k k k k k k k Kk k k K k k k k k k k k k k k k k k k k k k k k k k k K k k k k K k k k k k k k k k K k k k k k k k k k k k k k k kO h k KO KO KO RO K k h k h K kh K k x 


CALL THE APPROPRIATE EVENT SUBROUTINE 
0 0د ا‎ IASI IAAI IAI I IAAI IAAI IAI IAA IIA AIA IAAI IA ASAI IAAI AAAS ACAI A N 
40 CALL BRKDN 
GO TO 30 
50 CALL REPAIR 


Ya Wc Wr d dee ИИ Te de de de de Te We de de de e Te dae Te de r ohr kr r r We de de de de Te r de de Te e r e de de or ir r r r r kr de d e d kr de de r r r r r kr r r r r r r r r r r r r or r r kr r r r or r kr ok SSS. 


GO TO TIME ADVANCE ROUTINE TO SEE IF SIMULATION IS COMPLETE 


TT k e k k ft k k k r kr rr k k k k k k K k kk K K K K K K K h k K K k k k K k K kk K K k r or k k k k k k k k k kt k k k k k kt k k k r de k k k k k k k k K kt k k k k K k K K de k r ok or o o e e ŵr ŵr ŵr ŵr fn 


GO TO 30 


Te fe fe Te d d Te de krik k k k k N K K k ft K k K de kk kt de k k de de k k K K k k kt k kt k k k k k k k k kt k kt kt k k k k k k k k k k k k k k k k k k kt k k k k k k kt kt k k k k k k KK K k k k k K k k k k k k k k n Hœ 


WHEN SIMULATION IS OVER CALL REPORT GENERATOR 
KAK KKK KAKA KA NAA NA KAK NAA NN KAA AKKA KAKA KN KA KKK SC N KAA KAKA KA KAN KKK KA N NK OK KOK KOK K KOK KOK KOK He de Re KOR 
60 CALL RPTGEN 

STOP 

END 


SUBROUTINE INITLZ 


fe fe fe k K K K kk k k KK de De k k k k K K KK K KK k K K K k KO K k k K k k KO de de de de de de de ŵr de de de r or ok or or ok ok or r ok or or ok or r oir r or vir or or or r r K k K k k k KO ŵr rr k k k k h k k k k ŵr k ŵr ŵr ŵr Ŵr ŵr ŵr ŵr ŵ ŵ Xx 


INITIALIZATION ROUTINE 
LE E E E 5 8 5 5 55 E 55 5 555 585 8 8 5 E 5 E 55 E 55 5 8 8 5 5 8 5 5 8 8 5 5 8 5 85 5 5 5 5 5 8 5 5 8 3 5 24 3 8 5 85 5 85 5 5 Z 8 8 E 3 E 5 4 Z 5 5 8 EEL 35 25 24 5 45 2424385424452542 K k k n r 
INTEGER F,ELAWMT(100),TNREP 
REAL MITOF1,MITOF2,MITWS1,MITWS2,MITPT1,MITPT2,MITPC1,MITPC2, 
1 MSTOF1,MSTOF2,MSTWS1,MSTWS2,MSTPT1,MSTPT2,MSTPC1,MSTPC2, 
2  FEL(17),CHKIN(8), CHKINQ(100),NORML , TREPTM(B) , TTMINQ(B) 
COMMON /SIM/ MITOF1,MITOF2,MITWS1,MITWS2,MITPT1,MITPT2,MITPCL, 
1 MITPC2,MSTOF1,MSTOF2,MSTWSL,MSTWS2, MSTPT1,MSTPT2, MSTPCT,MSTPC2, 
2  SGMOF1,SGMOF2, SGMWS1, SGMWS2 , SGMPT1, SGMPT2,SGMPC1, SGMPC2,NAWTMT, 
3 NMTBUS, CHKIN,CHKINQ, TREPTM,MQ,F, TNREP, NMTPER,ELAWMT,TTMINQ, 
COMMON/TIMEKP/ CLOCK, IMEVT,NUMEVS, FEL 
Yr k k K k k k ir k k K k k K k k K k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k kt k k k f k k k k k k k k Kk k k k 
INITIALIZE SIMULATION 
1) SET SIMULATION CLOCK ТО 0 
2) ASSUME SYSTEM IS EMPTY AND IDLE AT TIME O 
3) INITIALIZE CUMULATIVE STATISTICS TO O 
* k'ik k K k kk k Kt k kt k k k de kk De kt k k k KK k KK k k k h h K k k k k k k k k k h k k k k k k k K k Kk k k k k k k k k ŵr de We de de k k K k k k k k k k k k k k k K k t k de de de de k k k k k K K kO k O k K k Wœ 
DO 10 I = 1,8 
TTMINOCI) 
TREPTM(I) 
10 CONTINUE 
CLOCK = 0.0 
IMEVT = 0 
NAWTMT 
NMTBUS 
TNREP = 0 
F=0 
MQ = 0 
fk k kk k K k k k k kk kk De k kt k k k k k k k k k k k k k k k k k k K de k k k k k k k de k k k k k k k k k k k k k k k k k k K e de de de de de de de kk de de de de de de de k K k k k de de K KO KO KOR ŵr ŵr ŵr ŵr ŵr Ŵr r ŵr n Hœ 
GENERATE TIMES FOR THE FIRST BREAKDOWNS AND SCHEDULE THEM IN 
FEL(1) - FEL(8), SET FEL(9) - FEL(16) TO INFINITY TO INDICATE 
THAT REPAIR IS NOT POSSIBLE WHILE THE SYSTEM IS UP, SET 
FEL(17) EQUAL TO LENGTH OF TIME SIMULATION IS TO RUN. 


Wrth k k ŵr k k de De De de e De De de de de de de de de fr cr rr dr cr dr de de de De de e de o r or kr o o ok r or o o ok ok o r r o kr r o e de de k k k k k K k k k k k k kk KO kO k KO k K kO K k k k k k kO k de de e e e ŵr ŵr ŵr 


0 
0 


0 
0 


FEL(1) = CLOCK * EXPON(MITOF1))**(1/.84) 
FEL(2) = CLOCK + EXPON(MITOF2))**(1/.84) 
FEL(3) = CLOCK + EXPON(MITWS1))**(1/.7) 

FEL(4) = CLOCK + EXPON(MITWS2))**(1/.7) 

FEL(5) = CLOCK + EXPON(MITPT1))**(1/.77) 
FEL(6) = CLOCK * EXPON(MITPT2))**(1/.77) 
FEL(7) = CLOCK + EXPON(MITPC1))**(1/.73) 
FEL(8) = CLOCK + EXPON(MITPC2))**(1/.73) 
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FEL(9) = 1.0E+30 
FEL(10) = 1.0E+30 
FEL(11) = 1.0E+30 
FEL(12) = 1.0E*30 
FEL(13) = 1.0E+30 
FEL(14) = 1.0Е+30 
FELC15) = 1.0Е+30 
FEL(16) = 1.0Е+30 
FELC17) = 6000 
RETURN 

END 


SUBROUTINE TIMADV | 
de Ye se de hr he e Ve dir dèr Vr dr dr dèr rf er rre e eee eee enne e e er n e e e n n n n e n n ne n n nn dèr de de dèr de dèr dèr e e n nee e e ne e n n e e n n š š š 3 
TIME ADVANCE ROUTINE FINDS NEXT EVENT ON FUTURE EVENT LIST 
AND ADVANCES CLOCK 
wr hr vh e he e err rr rr r dèr dèr vèr dr vc nt ete dc c dc dc e dc vc vc vi vk or vr oir vi or o vr r or or or or vr vr o vr vr vir or vr or vir vr vk r vr vir vr or yr vr or or or or vr or or vr or vr ok vr dèr e š Š 3 
INTEGER F,ELAWMT(100),TNREP 
REAL MITOF1,MITOF2,MITWS1,MITWS2,MITPT1,MITPT2,MITPC1,MITPC2, 
1 MSTOF1,MSTOF2,MSTWS1,MSTWS2,MSTPT1,MSTPT2,MSTPC1,MSTPC2, 
2 FEL(C17),CHKINC8), CHKINQ(100),NORML, TREPTM(8), TTMINQ(8) 
COMMON /SIM/ MITOF1,MITOF2,MITWS1,MITWS2,MITPT1,MITPT2,MITPCL, 
1 MITPC2,MSTOF1,MSTOF2,MSTWSL,MSTWS2,MSTPT1,MSTPT2,MSTPC1,MSTPC2, 
2. SGMOF1,SGMOF2, SGMWS1, SGMWS2, SGMPT1, SGMPT2, SGMPC1,SGMPC2,NAWTMT, 
3  NMTBUS,CHKIN,CHKINQ, TREPTM,MQ,F,TNREP,NMTPER,ELAWMT, TTMINQ 
COMMON /TIMEKP/ CLOCK, IMEVT,NUMEVS,FEL 
FMIN = 1.0E+29 
IMEVT = 0 


t she t k k k k k t k k k Yr k k k k k k kk k kk k k fr k k k k kk k k r k k k k k k k kk k k k k k k k k k k k k k k k k k k k k k k k k k k k k kk k k k k k kk k k k k k k k k k k k k k k k r r k r k K k e e e 


SEARCH THE FUTURE EVENTS LIST FOR THE NEXT EVENT 
k k w k k k W W k k k W k W k k k k k k k k k k k k k k k k k k k k k k k k k k k k. K k k k k k lr k k k k k kk k K kk k k kk k k k k k k k k k k k k kk k k k k k k k k k k k k k k k k k K K k k k k k À À K 
DO 30 I=L,NUMEVS 
IF (FEL(I).GE.FMIN) GO TO 30 
ЕМІН - ҒЕ((1) 
ІМЕУТ = I 
30 CONTINUE 
IF (IMEVT.GT.0) GO TO 50 


ft she he de k W k k k f k k k k k k k k k k k k k k W k k k k k k K k K k k r W k k k k dèr k k k k k k k r k k K k k k k k k k k k k k k kk k kk K KKL LLK k k k k k É k É k k k de de de de e de de e e œ œ œ x 


ERROR CONDITION - FUTURE EVENTS LIST IS EMPTY 
k e fe W k k k kk k k k k k k r k k k k k k k W k k k k k k W r k k k k k k fr k k k k k k t k t k k k k k k t k k k k k k k k fr k k k k Ék k k k k k k Ék k k k k k k k k k k k t k k k k k k k k k Kt Kk x 
40 — FORMAT(2X, '***** FUTURE EVENTS LIST EMPTY*****', 
C/,1X,'***** SIMULATION CANNOT CONTINUE *****') 
CALL RPTGEN 
STOP 
We vk vk k k e ei e e e k k k k k Ék k Ék k k k k k k k k k Ék k c e n k k k k k n e rne c n ree n d e x 
ADVANCE SIMULATION CLOCK 
NEXT EVENT IS TYPE 'IMEVT' WHICH WILL OCCUR AT TIME FEL(IMEVT) 
* de de de de hn k k k k k k k k k k k k k k e k k k K k k k k k k k k kk k K k k k k k k k k k k k Kk k k k k k Ék k k k k k k k k k k k k k йа ЖҰТ e re n tr kr k e r k ex 
50 CLOCK = FEL(IMEVT) 
RETURN 
END 


SUBROUTINE BRKDN 
Wk k k k k k W k k k k k k k k k k fr W k W k fr k k Kk k k k k fr k k r k r k k Ék k k k k k k k k t k k k k k k dèr k k k k k k k k k Ék k k k k k k k k k k dèr dèr dèc dèr dèr k kr r kt kr k t e k r t de de k r t k r r 
BREAK-DOWN EVENT ROUTINE - SCHEDULES THE BREAK-DOWN EVENT 
MAINTENANCE, OR TO A QUEUE AWAITING MAINTENANCE, COMPUTE 
MAINTENANCE TIME AND UPDATE STATISTICS. 
We de de k k k kk k k k Yr i Wr k W k k r k W W K k K k W K K K K K W K K K k K k k K k Kk K K kk k k k k k k k Kk k k k k k k k k k k k k k k k k k r vr vr tr $r K k r t Kk r r t k n o o o 
INTEGER F,ELAWMT(100),TNREP 
REAL MITOF1,MITOF2,MITWS1,MITWS2,MITPT1,MITPT2,MITPC1,MITPC2, 
1 MSTOF1,MSTOF2,MSTWS1,MSTWS2,MSTPT1,MSTPT2,MSTPC1,MSTPC2, 
2 FEL(17),CHKIN(8), CHKINQ( 100) , NORML, TREPTM(8), TTMINQ(8) 
COMMON /SIM/ MITOF1,MITOF2,MITWS1,MITWS2,MITPT1,MITPT2,MITPCL, 
1 MITPC2,MSTOF1,MSTOF2,MSTWS|,MSTWS2,MSTPT1,MSTPT2,MSTPC1,MSTPC2, 
2 SGMOF1,SGMOF2, SGMWS1,SGMWS2, SGMPT1, SGMPT2, SGMPC1, SGMPC2, NAWTMT, 
3 NMTBUS,CHKIN,CHKINQ, TREPTM, MQ, F , TNREP , NMTPER, ELAWMT, TTMINQ 
COMMON /TIMEKP/ CLOCK, IMEVT, NUMEVS, FEL 
SVT = 0 


Me de k W k k k k k k k k k k k k k kr ir i k r k k k r k k e e k k k k k k k k k k e k r K k K k k k k k 0 HHH n d x 
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SCHEDULE APPROPRIATE TYPE OF BREAK-DOWN 

20 FOR OFT1 BREAK-DOWN 
30 FOR OFT2 BREAK-DOWN 
40 FOR WST1 BREAK-DOWN 
50 FOR WST2 BREAK-DOWN 
60 FOR PTT1 BREAK-DOWN 
70 FOR PTT2 BREAK-DOWN 
80 FOR PTC1 BREAK-DOWN 

90 FOR PTC2 BREAK-DOWN 


+k k k krik حزان‎ fr c Te de e De De Ne e He dc de He de Te de De de Ve de Ve De Ve Te de e de Te de De de de e ee e ee k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k K k k k k k k k k k k k k k k 


GO TO (20,30,40,50,60,70,80,90),IMEVT 


Yet k k k k k t k t kt k k k k k k t kt k k k k k k k k k k e ir fr k k k k K k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k Kk ŵr Wr жж Ж 


DETERMINE IF MAINTENANCE PERSONNEL ARE ALL BUSY 


yrr i ir ic ŵr ŵr ŵr ŵr ŵr КККК К ККК اد بأد حفر‎ r ŵr ŵr KK KK e de âr e cr âr de die de de de de de Te Te Te de de Te de Te e e de de kr ir ir kr rr dr ir die de de e rc irr ŵr die de de de de ŵr ir ŵr ŵr e e ŵr de ŵr e ok 


20 IF(NMTBUS.GE.NMTPER) GO TO 25 
He he ke ke e de rr r de de e Ae e dr de de de ede de de de de de ke ke keke kede dede rde he e e ke de de kr de de ke e e e e e k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k 
MAINTENANCE PERSONNEL ARE AVAILABLE, UPDATE SYSTEM STATE AND 
COMPUTE SERVICE TIME FOR REPAIR 
ЖЕКИ ККК КККК ИКИ КИТКИ ИИТИИ КИК ИИВ ИИИ КИИ И ИИИ И ИИИ КИ ИИИ КИКИ КУК КККК КИК ИКИ И ИККИ К ИИК ИИ КИКИ КИКИ И ir âr de de Ve kr Ve de Âr ŵr ŵr ŵr e ŵr 
NMTBUS = NMTBUS + 1 
CHKIN(l) = CLOCK 
SVT = NORML(MSTOF1, SGMOF1) 
FEL(9) = CLOCK + SVT**(1/.33) 
FEL(1) = 1.0E*30 
GO TO 26 
Tc ie КК ŵr ŵr ic ic rc cc rc ic rc cr cc cr cic ir de de Ve Ж йй йй e e De e de ir de e de de e e e e De de de e e e e e e e e De e e e kr de e de de dde e e die Ve Ve e ee te e e e de de de de de de e e e d e x 
MAINTENANCE PERSONNEL NOT AVAILABLE, PUT REPAIR IN THE Q 
UPDATE STATISTICS 
+ k k k k k k k k'ik k k k k k k e k k k t k k k k k k k k k k k k k k k K k k k k k k k k k k k k k k k k k k k k k k k k k K k k k k k k k k k k k K k k k k k k k k k k k k k k k k k k kt k k k k k k 
25 NAWTMT = NAWTMT + 1 
ع1‎ (NAWTMT.GT.MQ) MQ = NAWTMT 
IF (NAWTMT.GT.99) GO TO 100 
CHKINQ(NAWTMT) = CLOCK 
ELAWMT(NAWTMT) = 1 
FEL(1) = 1.0E+30 
26 RETURN 


Tc fe he e he e De Te Ve e de e e de De Ve Ve de e Ve Ve Ve Ve Be Ve Ve Ve Ve Ve Be Ve de Ve Ve De Ve e de Ve de de Ve de de e de de de de Ve de ir fr rr rr rc cr rr ŵr rir De de de De e de de Yr cr rrr rr Ve ŵr e de Ve de Ve Ve de de M ir de de Wœ 


DETERMINE IF MAINTENANCE PERSONNEL ARE ALL BUSY 


*fe e he he e hee e e de e Te de K KKK KKK RAK KKK KRK KKK K KKK e e e e e e e e e eee e e e ee ce e e eh e e e e e e e e e e e de e ee ee ce ee de ie Be die de e ir ŵr Âr ŵr ŵr Âr ŵr ŵr e 


30  IF(NMTBUS.GE.NMTPER) GO TO 35 


fe he e e e e e ККК e ee ee e e ee ee e Ve he ehe e e e De e kr ok kr ok ok kr kr e e ee eee e kr ok kr kr ok ok e e e kr ok ok r kr e Ve e e e e e e Ve e de e Ve de Ve e Ve Ve e ok ok kr e e e e e e de Me We e Me x 


MAINTENANCE PERSONNEL ARE AVAILABLE, UPDATE SYSTEM STATE AND 
COMPUTE SERVICE TIME FOR REPAIR 
rc ŵr ŵr rc e ee hl e k k k k k k k k k k r r kr ce ee t k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k x 
NMTBUS = NMTBUS + 1 
CHKIN(2) = CLOCK 
SVT = NORML(MSTOF2,SGMOF2) 
FEL(10) = CLOCK + EXP(SVT) 
FEL(2) = 1.0E+30 
GO TO 36 
** kk k k k k k k k ik kk k Ne k k k k k k k t k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k 
MAINTENANCE PERSONNEL NOT AVAILABLE, PUT REPAIR IN THE Q 
UPDATE STATISTICS 
* k k k * ee le k k k k k k k k k k k k k k k k k k k k k k k k k k k k Ŵr k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k ike k e k k k k k k k k ŵr k k ŵr k k k ke 
35 NAWTMT = NAWTMT + 1 
IF (NAWTMT.GT.MQ) MQ = NAWTMT 
IF (NAWTMT.GT.99) GO TO 100 
CHKINQ(NAWTMT) = CLOCK 
ELAWMT(NAWTMT) = 2 
FEL(2) = 1.0E*30 
36 RETURN 


* k k k k k'ik kt ir Yt r k k ke k k ir Tt r k e ir ir k k k k k k k k k k k k k k k k ik k k k k k k e k k k ie k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k ŵr rr ŵr ŵr ŵr ŵr ŵr 


DETERMINE IF MAINTENANCE PERSONNEL ARE ALL BUSY 


Tek k k k k Tk Se e kt kt k k e k k k k k k k k k k k'ik k k k k k k k k k k k k k k k k k k k k k k k k e k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k de c de k Ŵ 


40 IF(NMTBUS.GE.NMTPER) GO TO 45 


ккк k k k k k k k k Tk k k k kk k k e ke k k ke k k k K k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k K fe fe We x 


MAINTENANCE PERSONNEL ARE AVAILABLE, UPDATE SYSTEM STATE AND 
COMPUTE SERVICE TIME FOR REPAIR 


45 


Wr tr ŵr k K К e ehe K ŵr rr câr rr ŵr Ŵr rrr ŵr Kk k КК جف‎ ke k he k k k K K K K k K K K KK K k r o k ir rr k k K k K k k k K k KK K k K k k k k K K k K Kk k K k k k k K K ŵr ir ŵr k Кї Кї К К yr 


NMTBUS = NMTBUS + 1 
CHKIN(3) = CLOCK 
SVT = NORML(MSTWS1,SGMWS1) 
FEL(11) = CLOCK + SVT**(1/.32) 
FEL(3) = 1.0E+30 
GO TO 46 
hr hr tèc tèc KK KAKAK KK KKK Ñ KAKAK KAK KA KK KKK KKK KKK KK KK KKK KKK KKK KKK KKK KKK KKK KK KK KKK KAKAK AKKA KKK KKK K KK KK K 
MAINTENANCE PERSONNEL NOT AVAILABLE, PUT REPAIR IN THE O 
UPDATE STATISTICS 


rr ŵr ŵr rr ŵr ir ŵr ŵr Ve Ye De Ye ŵr e cr ŵr de e rrr rr rh cce e cr c rc e d ir cr câr rr cr ir cr ŵr câr cr âr rr ŵr ŵr rr rr ŵr cr rr ir de ŵr ir ŵr ŵr ŵr cr ŵr ir cr die de die de de de de de e e n d d d d x 


45. NAWTMT = NAWTMT + 1 
IF (NAWTMT.GT.MQ) MQ = NAWTMT 
IF (NAWTMT.GT.99) GO TO 100 
CHKINQ(NAWTMT) = CLOCK 
ELAWMT(NAWTMT) = 3 
FEL(3) = 1.0E+30 

46 RETURN 


dc e k k e e ee nr rr rh ce ccc e cc c ir ŵr ŵr ŵr ŵr K K k kk k k k K k ŵr ŵr ŵr k k Kk Kk Kk Kk Kk ir ŵr ŵr ŵr ir k K K KK KK KK k kk k k K kk k k k k k k Kk k k r r ŵr ŵr ŵr ŵr ŵr câr cr cr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr 


DETERMINE IF MAINTENANCE PERSONNEL ARE ALL BUSY 


fk kk k k kk k K k k kt k kt k K kk k k k k k k Kk K K K kk k kk kk k k k kk k k k k k ŵr k ŵr ŵr k ŵr k k KK k Kk kk k k k ŵr ŵr ŵr ŵr ir k K k k k k k k k k k k k k k k ŵr k K Kk k k k k k Kk K K K K k k k k K k ŵr k K K 


50 IF(NMTBUS.GE.NMTPER) GO TO 55 
F A AAA A A A A A AA AA A A A A AAA AAA A AA AA AAA A fcre ef rr rre cfr rn rrr hr n rr e xn YK 
MAINTENANCE PERSONNEL ARE AVAILABLE, UPDATE SYSTEM STATE AND 
COMPUTE SERVICE TIME FOR REPAIR 
"dc Tc Uc Tc Tc e c de Se e Dc de de de de de de c îr ŵ ir ŵ * 11111 жая 
NMTBUS = NMTBUS + 1 
CHKIN(4) = CLOCK 
SVT = NORML(MSTWS2,SGMWS2) 
FEL(12) = CLOCK + SVT**(1/.32) 
FEL(4) = 1.0E+30 
GO TO 56 
*dr e e ЕЖ ЖЕ Х.Х Жи ж ж жж e d xn x 
MAINTENANCE PERSONNEL NOT AVAILABLE, PUT REPAIR IN THE O 
UPDATE STATISTICS 
xc d Tc Te Tc Te de de ЕКЕ de de Kk Kk de k K Kk Kk de Kt de de de Kf de de de de Kt Kk Kk f k k kt k k k Kk Kk K kt Kk k k K k kt Kk Kk Kk kt K Kt Me tf tf K K K K k k K k Kk Kk k Kt Kt Kt kf k Kk k Kk Kk k K kt kk Kk K K£ de Xœ K K Kk Kk K Wœ 
55 NAWTMT = NAWTMT + 1 
IF (NAWTMT.GT.MO) MO = NAWTMT 
IF (NAWTMT.GT.99) GO TO 100 
CHKINQ(NAWTMT) = CLOCK 
ELAWMT(NAWTMT) = 4 
FEL(4) = 1.0E+30 
56 RETURN 


e e e ŵr Wr c de de de Кї КККК de de ŵr de de ŵr de O de e de de de de e c de e rr ŵr rr ir cr cr ir ir câr ir ir câr câr cr cr cr cr rr cr cr ŵr die cr ŵr rr ŵr ŵr e de e die ŵr rr rr Ne e de he de Re Ne Me Me Me Ne Me N Me N e de de de ŵr ŵr d d e e 


DETERMINE IF MAINTENANCE PERSONNEL ARE ALL BUSY 


* fe k k K k k fk k k k k K k k t k k K k K k K k k k k kk ŵr ŵr ŵr K k kk Kk de de de de de e ee e e c ree e k ir ŵr ŵr cr ŵr ir cr rr ŵr ir ŵr rr rr ŵr rr ŵr rr K K KK K k k kk k k k k KK K KK K k k k k k k r e 0 ŵr ir À 


60 — IF(NMTBUS.GE.NMTPER) GO TO 65 
KKK KK ANAKE e De e De KAKAK KAK KAKAK KKK KKK KKK KAKA KK KKK KK KAAKA KAKAK Yk II I i ŵr ŵr ŵr ŵr rr e 
MAINTENANCE PERSONNEL ARE AVAILABLE, UPDATE SYSTEM STATE AND 
COMPUTE SERVICE TIME FOR REPAIR 
9 I IE ICT TIC De e de De De e de de di iii KAAKA AAA KAK KKK AA KA KKK KAKAK KAAKA A KKK KKK KKK KAKAK KKK KKK KKK KA KAK KKK KK KAKAK KA 
NMTBUS - NMTBUS * 1 
CHKIN(5) = CLOCK 
SVT = NORML(MSTPT1,SGMPT1) 
FEL(13) = CLOCK + SVT**(1/.5) 
FEL(5) = 1.0E+30 
GO TO 66 
KAKA KKK KAAKA AKKA KAA KKK eee eee dede del do dede KKK KKK A KKK KAKAK KAKAK KAKAK KKK e de de de de Te de dh de de de nit tt ttt 
MAINTENANCE PERSONNEL NOT AVAILABLE, PUT REPAIR IN THE O 
UPDATE STATISTICS 
KKK A AAA AKA K HR RR RR RR rrr rrr ddr rrr rrr rrr dde ooi rr rrr rr ddr rrr rr rr rrr rr ŵr rr 
65  NAWTMT = NAWTMT + 1 
IF (NAWTMT.GT.MQ) MQ = NAWTMT 
IF (NAWTMT.GT.99) GO TO 100 
СНКІМО(МАЫТМТ) = CLOCK 
ELAWMT(NAWTMT) = 5 
FEL(5) = 1.0E+30 
66 RETURN 
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Pr tc e e de ŵr ŷr ŵr ŵr ŵr ŵr ŵr ŷr ŵr ŵr de de de de De de de de de ce de k de de de de de de de de de ŵ * de de de de de kc de Pe de de de de Me de de de * ŵ * ŵ * ŵ * ŵ ŵ ŵ * de ŵ ŵr H e ŵ ŵ ŵ * ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵr ŵr  ŵ ŵ ŵ ŵ ŵ * * r 


DETERMINE IF MAINTENANCE PERSONNEL ARE ALL BUSY 


Te fe fe e dede de de dee dee De De De Oe de de e de КККК ККК И Ŵr Ŵr Ŵr ŷr dr Ŵr de ŵr ŵr ŷr ŵr ŵr ŵr ŵr ŵr ŵr ic rr de de de de de de de D de de de de fc de de e ŵr ŵr dr fr ŵr dr ŵr de de de de de de de de de de de de de de ŵr rr de i ok o ok ok ok Yr yr r o o ok 


70 IFCNMTBUS.GE.NMTPER) GO TO 75 
re 0 d d Se UR WR HB RH HR جا جا‎ eee dede De de de de de de eee He He He He de de de de de de de RR RW de de de de de d de de de H de de de de D H KAKAK AKKA KK A 
MAINTENANCE PERSONNEL ARE AVAILABLE, UPDATE SYSTEM STATE AND 
COMPUTE SERVICE TIME FOR REPAIR 
> he 9 3% dc Se de De De hr e De Te Te Te De e de OR lode de de: de de de ITT oe eec doc dedo Ne de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de KO KORR R R R R œ 
NMTBUS = NMTBUS + 1 
CHKIN(6) = CLOCK 
SVT = NORML(MSTPT2,SGMPT2) 
FEL(14) = CLOCK + SVT**(1/.5) 
FEL(6) = 1.0E+30 
GO TO 76 
KAAKA AAA KA AAA KA KK AKKA AAA KKK de Ne de Ne Hd He de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de de e % e Ke 
MAINTENANCE PERSONNEL NOT AVAILABLE, PUT REPAIR IN THE Q 
UPDATE STATISTICS 
AAA AKA AAA A KAA KAKA KKK KK AAA A KANA de De de de de de ink 
7959 МАМТМТ = МАМТМТ + 1 
IF (NAWTMT.GT.MO) MQ = NAWTMT 
IF (NAWTMT.GT.99) GO TO 100 
CHKINQ(NAWTMT) = CLOCK 
ELAWMT(NAWTMT) = 6 
FEL(6) = 1.0E+30 
76 RETURN 


fer e e he dee ede ede de de de KR k Me k de k k k k k k k k k k k k Kk k k k kt k k k k k k k k kk k k k k k k k k k k kk k k k k k k k k k kk k k k k k k k k k kk k k k k k k k k k kt k k k k k k k k k k Wœ 


DETERMINE IF MAINTENANCE PERSONNEL ARE ALL BUSY 


de I de de de de e ŷr cr ŵr Ue de de de e de de e de de ŵr de Me de de Me r o o kk i vk r r r YY ok oir r Yk ok ok o Yk  ŵ ŵ ŵ ŵ ŵ ŵ Me ŵ ŵ ŵ ŵ ŵ e de ŵ ŵ ŵ de ŵ ŵ ŵ ŵ ŵ ŵ ŷr ŵ ŵ ir ŵ ŷ de ŵ ŵ ŵ de de de de de de H e ir ok r yr oir ok oir ok o ir ok 


80 — IF(NMTBUS.GE.NMTPER) GO TO 85 
VA AAA AA A A A AA AA A A Me A A e جلا جا‎ de de ŵ de ŵ de de Me de Me ŵ ŵ  ŵ ŵ  ŵ 'k ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ De Me de ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ir ŵ ŵ ŵ ŵ  *  ŵ ŵ ŵ y ŵ ŵ i ŵ * ŵ ŵ 
MAINTENANCE PERSONNEL ARE AVAILABLE, UPDATE SYSTEM STATE 
COMPUTE SERVICE TIME FOR REPAIR 
kr fede dede de e de De de Ie de de de de de de de de de de ج11‎ de de de de de de e de de de de de e de de de de de de Me de de de de de de de de de de de de de de de e de de de de de de Me de de Me de We de de e de de de de de de de de de de de de de de r r i i r Yr yr yr yr yi o r oi 
NMTBUS = NMTBUS + 1 
CHKIN(7) = CLOCK 
SVT = NORML(MSTPC1,SGMPC1) 
FEL(15) = CLOCK + SVT**(1/.14) 
FEL(7) = 1.0E+30 
GO TO 86 
La AAA AAA AA AAA A de e de de de de de de de de de de de de de de e ج11‎ de de de de de de de de de e de e e de de de de de Me de de de de de de de de de e de de de de de de de de de de de de de de de de e eee eee r vr ir oi yr oi r oir yr o d x 
MAINTENANCE PERSONNEL NOT AVAILABLE, PUT REPAIR IN THE O 
UPDATE STATISTICS 
LA AA AAA de de جلا‎ e e e جل‎ de de de De Pe de جلا جلا‎ de  ŵ ŷ * ŷ de ŵ ŵ ŵ ŵ Me ŵ ŵ ŵ جلا‎ ŵ ŷ ŷ ŵ ŵ ŵ ŵ ŷ ŵ ŵ  ŵ ŷ de Me de ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŷ i de ŵ ŵ ŵ ŷ r de De ŵ de  ê ŷ * ŷŵ ŵ ŵ * ŵ ŵ  *  ŵ * * * ŵ ŵ ŵ'k e e d n x x 
85 NAWTMT = NAWIMT + 1 
IF (NAWTMT.GT.MO) MO = NAWTMT 
IF (NAWTMT.GT.99) GO TO 100 
CHKINQ(NAWTMT) = CLOCK 
ELAWMT(NAWTMT) = 7 
FEL(7) = 1.0E+30 
86 RETURN 


125444444424 GG k k k k k k k k K k k k k k k kt k k k k k k k k k k k k K K K k k k k k K Kk K k k k k k k k k k k k k k k K k k k k k k k k k k k k k kt k k k k k k h.h 


DETERMINE IF MAINTENANCE PERSONNEL ARE ALL BUSY 


de de de dede de de e de e Ie Ye de e Ïe De e de e de de de De e de de de de de de de de de de de e de De de de de Me de جا‎ de de de de de de de de de de de de Kh de de de de de de de de de de de e de de de de de de de de de de de de H de R He R R Œ H H e ee ded dex 


90 — IFCNMTBUS.GE.NMTPER) GO TO 95 
de Se dc e e e % IIIT IIT ŵr rr rr ŵr k K K K K K K de k K k K K K K K K KOK K K K K K K K K de de de de de de de K de de KOK K K K K K K K de de de K K K K de de de de de de de de de K K K 
MAINTENANCE PERSONNEL ARE AVAILABLE, UPDATE SYSTEM STATE AND 
COMPUTE SERVICE TIME FOR REPAIR 
de de dc 3% He KAAKA UU UD de de de de de de de de de de d de d de de de de de de eee eode de deed eee dee КЕКИ К 
NMTBUS = NMTBUS + 1 
CHKIN(8) = CLOCK 
SVT = NORML(MSTPC2,SGMPC2) 
FEL(16) = CLOCK + SVT**(1/.14) 
FEL(8) = 1.0E+30 
GO TO 96 
ÁR SCC III IO OOO de de De de de De H Se de de جلا‎ de de de de de de de de d de de de de de d d He d d dede dede de He d de de de de d de de de de de de de de de d de de de de d de d d R R d eek 
MAINTENANCE PERSONNEL NOT AVAILABLE, PUT REPAIR IN THE Q 
UPDATE STATISTICS 


* kk k k kk kt k k k kt k k k k k k k k k k k k k k k k k kk k k k k k k k k k k k k k k k kk k k h k k k k k k k k k k k K k kk k k k K k K K h R Kk k de K k Kk k K Kk k K AA K 
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95 NAWTMT - NAWTMT + 1 
IF (NAWTMT.GT.MO) MO - NAWTMT 
IF (NAWTMT.GT.99) GO TO 100 
CHKINQ(NAWTMT) = CLOCK 
ELAWMT(NAWTMT) - 8 
FEL(8) = 1.0E+30 
96 RETURN 
100 WRITE(6,105) 
105 FORMAT(3X,'***** TO MANY AWAITING MAINTENANCE ACTIONS *****! /, 
C1X,'***** INCREASE NUMBER OF MAINTENANCE PERSONNEL *****') 
CALL RPTGEN 
STOP 
END 


SUBROUTINE REPAIR 
ŵr îr ŵr ŵr hr КЛ Ŵr ŷr Ŵr ŷr ŷr ŷr ŷr Ŵr ŷr Ŵr ŷr Ŵr Ŵr ŷr ŵ ŷr Ŵr Ŵ ŷr ŷr Ŵr ŷr Ŵ ŵr ŷr Ŷr ŵ ŵr ŷr ŵr ŵ Ŵr ŵr ŷr ŷr Ŵr Ŵ ŷr Ŵr îr Ŵr ŵr ŷr Ŵr Ŵr ŷr ŵ ŷr ŷr ŷr Yr Ŵr ŵ ŵ ŷr Ŵr ŷr Yr Ŵr Yr Ŵ ŵ ŷr ŵr ŵr ŷr ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ n xn x 
REPAIR EVENT ROUTINE: NOTE REPAIR COMPLETE, TAKE NEXT 
BREAK-DOWN EVENT OUT OF THE QUEUE, COMPUTE REPAIR TIME, 
SCHEDULE REPAIR EVENT. 
ŵr ŵr ККК ŷr ŷr ŵr ŷr ŷr ŵr ŷr Ŷ ŵ ŷr ŷr ŷr ŵ ŵ ŵ ŷr ŵ ŷr ŵ ŵ ŷr ŵ ŷr ŵ ŵ ŵ ŷr ŵ ŷr ŵ ŵr ŵ ŷr ŷr ŵ ŷr ŵr ŵ ŷ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŷŵ ŵ ŵ ŷŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŷr ŵ ŷr ŵ ŷr ŵr ŷr ŷr ŵ ŷr ŵr ŵ ŵ ŵ ŵ ŵ ŵ ŵ 
INTEGER F,ELAWMT(100),TNREP 
REAL MITOF1,MITOF2,MITWS1,MITWS2,MITPT1,MITPT2,MITPC1,MITPC2, 
1 MSTOF1,MSTOF2,MSTWS1,MSTWS2,MSTPT1,MSTPT2,MSTPC1,MSTPC2, 
2 FEL(17),CHKIN(8),CHKINO(100),NORML,TREPTM(8),TTMINO(8) 
COMMON /SIM/ MITOF1,MITOF2,MITWS1,MITWS2,MITPT1,MITPT2,MITPCL, 
1. MITPC2,MSTOF1,MSTOF2,MSTWSL,MSTWS2,MSTPT1,MSTPT2,MSTPC1,MSTPC2, 
2 SGMOF1,SGMOF2,SGMWS1,SGMWS2,SGMPT1,SGMPT2,SGMPC1,SGMPC2,NAWTMT, 
3 NMTBUS,CHKIN,CHKINQ, TREPTM,MQ,F, TNREP,NMTPER, ELAWMT, TTMINQ 
COMMON /TIMEKP/ CLOCK, IMEVT,NUMEVS, FEL 
de de ŷr ŵr ŷr ŵr ŷr ŵr ŵr ŷr Ŵr ŵr ŵr ŷr ŷr ŵr ŷr yrsŵr Ŵr ŵr ŵr ŷr ŷr ŷr Ŵr ŷr ŷr ŷr Ŵr ŵr ŵ ŵ ŷr ŷr ŷr ŵr ŵr ŷr ŷr ŷr ŷr a o 
SCHEDULE APPROPRIATE TYPE OF REPAIR 
20 - OFT1 REPAIR 
30 - OFT2 REPAIR 
40 - WST1 REPAIR 
50 - WST2 REPAIR 
60 - PTT1 REPAIR 
70 - PTT2 REPAIR 
80 - PTC1 REPAIR 
90 - PTC2 REPAIR 


ŷr ŷr ŷr ŷr ŷr ŷr Ŵr De ŵr ŵr ŵr ŷr ŵr yr Me ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŷr ŵr ŵr ŵr ŵr ŵr ŷr ŵr ŵ ŷr ŵr ŵ ŷr ŵr ŷr ŷr ŷr ŵ ŵr ŵr ŵr ŵr ŵr ŵr ŷr ŵ ŷr ŵr ŷr ŵ ŷr ŵr ŵ ŵ ŵ ŷr ŵ ŷr ŷr ŷr ŵr ŷr Yr ŵr ŵ ŵr ŷr ŵ ŵr ŷr ŵ ŷ ŵ ŵr ŵ ŵ ŷr ŷr ŵr ŷr ŵr ŷr ŷr ŵ ŵ ŷr ŷr ŵ ŷr ŵ ŵ ŵ ŵ ŷr ŵ ŷr He 


GO TO (20,30,40,50,60,70,80,90), IMEVT-8 


Yr Yr ŵr ŵr ŵr ŷr ŷr ŵr ŵr ŷr ŷr ŵr ŷr ŷr ŷr ŵr yr ŵr ŵr ŵr ŷr ŵr ŵr ŵr ŵ ŵr ŷr ŷr ŵr ŷr ŵr ŷr ŵ ŷr ŵr ŵ ŷr ŵr ŵ ŵ ŷr ŷr ŵr ŵr ŵ ŵr ŵr ŷr ŷr ŵ ŵ ŷr ŷr ŵr ŵ ŷr ŷr ŵ ŵ ŷr ŵ ŵ ŷr kk k ŵr ŵ ŷr ŵ ŷr ŵr ŵ ŷr ŵr ŵ ŵ ŵr ŵ ŵr ŵ ŵ ŵr ŵr ŷr ŷr ŵ ŵr 


UPDATE CUMULATIVE STATISTICS 
Ye scd de de jejeje je de de de dee dee de c de de dee de ee de e fe de e fe dee de de e e dece de e c de e de de e de e e de ee e ee de dee de e e dee ede e je e e o ŵr rr e 
20 TREPTM(1) = TREPTM(1) + (CLOCK - CHKIN(L)) 
TNREP = TNREP + 1 
IF((CLOCK - CHKIN(1)).GE.3.0) F = F + 1 


dr ŵr yr Ŵr ir k i ir k e ŵr ŷr ŷr ŵr ŷr Ŵr t Wr Ŵr rr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŷr ŵr ŵr ŵr ŷr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr Ŵr ŵr Ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr Ŵr Ŵr Ŵr ŵr ŵr ŵr ŵr ŵr ŷr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵ ŵr ŵr ŵr ŵr ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵr ŵ ŵr ŵ ŷr ŷr ŵ 


CHECK CONDITION OF OUEUE 


lt Se ŷr ŵr Ŵr Ŵr ŷr Ŵr rc ŵr yr sêr ŷr yr ŷr ŷr ŵr ŷr ŷr ŷr ŷr ŷr ŷr ŷr Ŵr ŷr ŷr Ŵr ŷr ŷr Ŵr ŷr Ŵr ŷr ŷr ŷr ŷr ŷr ŷr ŷr ŷr ŷr ŷr Ŵr e fe ŷr fe ŷr fe Ŵr Ŵr ŷr ŷr ŷr ŷr ŵr ŵr Ŵr ŷr ŷr ŵr ŵr yir ŵr Ŵr Ŵr ŵr ŷr ŵr ŵr ŵr ŷr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŷr kkk k ا ا‎ 


IF(NAWTMT).GE.1.0) GO TO 25 
hr Ŵr ŵr n e ŷr yr ŵr ŵr ŷr Ŵr sîr جارد‎ ŵr îr جاو‎ ŷr ŷr ŵr ŷr Ŵr ŷr ŵr ŷr ŵr ŷr ŷr ŷr ŵr ŷr ŵr Ŵr ŵr ŷr ŵr Yr ŷr ŷr ŵr ŷr ŵr ŷr Ŵr Ŵr ŵr ŷr ŵr ŷr ŷr ŵr ŷr ŷr Ŵr ŵr ŷr ŵr ŷr iir ŵr ŵr yir ŵr ŵr ŵr ŵr ŵ ŵ ŵr ŷr ŵr ŷr ŵr ŷr ŷr ŷr ŵr ŵr ŷr ŵr ŷr ŷr ŷr ŵr ŷr ŵr ŷr ŵ ŷr ŵ ŵ ŵ ŵr ŵ ŵ ŷŵ ŵ x 
QUEUE EMPTY, SERVICE MAN BECOMES IDLE, NEXT DEPARTURE TIME 
EQUAL TO INFINITY 
Ŵr ŵr ŵr ŵr de rrr rr sîr ŵr dr Ŵr ŵr ŷr ŷr rir ŵr ŷr ŵr ŵr جاو‎ ŵr ŷr ŵr ŵr ŵr ŷr ŷr ŵr ŵr ŷr ŵr Ŵr ŷr ŷr ŵr ŷr ir Ŵr Ŵr ŷr ŵr yr ŷr Ŵr ŵr ŷr ŵr ŵr ŷr ŵr ŵr ŵr ŷr ŷr ŷr ŵr ŵr ŵr ŵr ŵ ŵ ŵr ŷr ŵr ŵ ŷr ŵr ŵr ŵ ŷr ŵ ŵ ŷr ŵr ŵr ŵ ŷr ŵ ŵ ŷr ŵ ŵ ŷr ŵr ŵ ŵ ŵr ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵr ŵ 
NMTBUS = NMTBUS - 1 
FEL(9) = 1.0E+30 
GO TO 28 
ŵ ŷr ŷr ŷr ŵr ŵr ŵr ŵr ir ŵ ir fc ŵr ŵr ŵr Yr rir ŷr ŷr ŷr ŵr ŵr ŵ ŷr ŷr ŵr ŵr ŵr ŵr ŵ ŵ ŷr ŷr ŷr ŵ ŷr ŵ ŵ ŷr ŵ ŷr ŵ ŵ ŵ ŵ ŷr ŵr ŷr ŷ ŷr ŵ ŷr ŵ ŵ ŵ ŷr ŵr ŵ ŷr ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŷr ŷr ŵr ŵr ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵr ŵ ŵ ŵ ŵ ŷr ŷr Ŵr 
QUEUE CONTAINS AT LEAST ONE MORE MAINTENANCE ACTION, START WORK ON 
THE NEXT MAINTENANCE ACTION IN THE O, AND MOVE-THE REST UP ONE 
POSITION IN THE STACK 
ŵr ŵr ŵr ŵr ŷr Ŵr irr ŵr yrr ir ŵr Ŵr ŷr ŵr ŵr ir 1 Yr ŷr ŷr ŷr ŷr ŷr ŷr Ŵr Ŵr ŷr ŷr ŷr Ŵr ŷr ŷr ŷr ŷr ŷr Ŵr ŵr ŷr ŵr ŷr ŷr ŵr ŵr ŵr ŵr ŵr ŷr ŷr Yr ŷr ŷr ŵr ŵr Yr ŷr ŷr ŷr ŷr ŷr ŵr ŵ ŷr Ŵr ŵr ŷr ŵ ŷr ŷr ŵ ŵr ŷr ŵ ŷr ŷr ŵr ŷr ŷr ŵ ŵ Ŵr ŷr ŵ ŵ ŵ ŵ ŵ ŷr ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŵ ŷ ŵ ŵ 
25 GO TO (26,36,46,56,66,76,86,96),ELAWMT(L) 
26 CHKIN(1) z CHKINO(1) 
TTMINQ(1) = TTMINQ(1) + (CLOCK - CHKINQ(L)) 
DO 27 I = 1,NAWTMT - 1 
CHKINO(I) z CHKINO(I*l) 
ELAWMT(I) z ELAWMT(I*l) 
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27 CONTINUE 


ir hk r 4 $ 2 2 4.2 2.6.6 0.8 25.4 2.4 4 2.4 4 E 2.24 2 2 4 2 2 4 à 8.2 4 4. 5 8 8 8 r ir ir ir r i ir r kr kr r 8 8 5 8 8 5 5 r kr ir E he ir r kr r ir kr r he ee e e t ir e e e e ir e e ir ir e ir ir r EEE EE EE 


UPDATE LENGTH OF QUEUE 
ir iriri k k Sir ir ir ie hr ir fir kr dir fir sir k YK hk T de de de dee de de de k k k e dee de جف جلا‎ De de e k جلا‎ dr جلا جلا‎ e de جلا جف جف‎ e de e k De جلا‎ k e de e e جلا جف‎ de de de e e k k k k de جلا‎ de de Se De جلا‎ e e He Se جلا‎ dr dr dr dr e dr á e de de de e 
NAWTMT - NAWTMT - 1 
SVT = NORML(MSTOF1,SGMOF 1) 
ҒЕ((9) = CLOCK + SVT**(Ll/.33) 
RETURN 


ЖЕТЕ De ОУС de De de die De de de de du De De De de Dde Ye de de de de Óe De De De Ye De Ye de Ye De De Ye r Ye De r De r r r Ne r Ye De r de r de r De r de Ne De r ir Ne de e De de De de Se de e rá de Me de ir rá de Me árd irí 


SCHEDULE NEXT BREAK-DOWN FOR OFT] 
de dodet AAA AAA KAE KK KKK KAK AAA AA KKK KKK AK AAA K KKK KKK KAAKA KKK KKK AKKA KKK AKA KKK dede dedo dede de dede ee ei dee ede 
28 IAT » EXPON(MITOF1) 
FEL(1) = CLOCK + 1АТ%%(1/.84) 
RETURN 


тай ŵr r A a a k T îr ŷr fr k k k k k k k k k Kk Ŵr dr k k k ŵr k k kt de ŵr k k Ŵr dr k k k k k k Kk Kk k k k k k k k k k k k k Kk Ŵr rc Kk ŵr ic k k k kt kt Kk kt kt k k k k k k k k ŷr k k k k kt k kt k k k k k ŷr ŷr ŷr 


UPDATE CUMULATIVE STATISTICS 
حفر جد حار جاه عاد جار‎ R R A A A R TA e A a a A L fr de de de de Ŵr dc ŷr r ŵr dc r dc dc ic ŵr r ŵr Ŵr dc r fr br ŵr îr ŵr ŵr dc îr fc br îr îr br fc br ŷr îr br ŵr br br ôr ŷr îr br îr îr bc îr îr bc de br De fc îr îr îr Óe fc de de fc * De îr ŵr c * de de de de de de ŷr de fc r 
30  TREPTM(2) = TREPTM(2) + (CLOCK - CHKIN(2)) 

TNREP = TNREP + 1 

IF(CCLOCK - CHKIN(2)).GE.3.0) F =F + 1 


hr hr vk fr k kt k rk k k k ŷr ŷr k k k k k k k k k k k k k K k k Ŵr k k k k k kt dc dr chr ŵr ft kk kt dr Kt k k kt k k kt k kt k k k kt k k k k kt br k kt kt k k kt k k dr k kt k k t kt k kt de br dc ir fc ŷr k k ŵr de ŷr i ŵr 


CHECK CONDITION OF QUEUE 


ft fe de de de e kk K k k kk k K k ŵr k k k k k k k k Ŵr k k k kt fr k kt k k k k kt k k k k k kr kr k k k dr cr k k k k k kr k kr k k kk kk k k Kk k k k k kk k k k k k k k k k k kr îr k k k k k k k k h k h h K K EEE E 


IF(NAWTMT.GE.1.0) GO TO 35 
de de de de dede de de de ŵr ár ár á Ŵr dc dc fc dr ár á dc dc ŵr cr ic ir îr îr ic ir De dc dc Ŵr ŷr dr ŵr ir fc ŵr Ŵr Ŵr dr ŷr dr ŷr ŵr ir ir Ŵr ŷr ŵr ŷr ŷr ŷr ic ŷr îr ŵr ŷr ŷr ŵr ŷr ŷr ŵr ŵr ŷr ir vir fr ŵr ŵr ŷr dr ŵr ŵr ôr ŷr ŵr r dr ŷr rr ir dr * ŵr ir ŵr ŷr ŵr ŵr ŵr ŵr ŵr x 
OUEUE EMPTY, SERVICE MAN BECOMES IDLE, NEXT DEPARTURE TIME 
EQUAL TO INFINITY 
de de de îr dc ic ŷr fc ŷr ôr de De br Se k r kt КЕ de k k k k k ŷr k k Kk kt kt k k k k ŷr k k k Kk k k k ŷr ŷr kt k k k k k k k k k kt k k kt k k k k k k ŷr ŵr kk k k k k k k k ir k kt Kk kt ŷr dc ft fc k ŷr ŷr k kt k ŵr ŷr ŵr ŵr Kk 
NMTBUS = NMTBUS - 1 — FEL(10) = 1.0E+30 
GO TO 38 


ŵr ŵr ŵr k k k k k k k kt k k k k k k kt k k k k k k k k k k k k kt k k k k k kt k k k k k k k k k k k k k k k k k k k k kt k k k kt k k k k k k k kt kt k cr k k k k k k k k k k k k k k k k k k k k k k k k f k k K 


QUEUE CONTAINS AT LEAST ONE MORE MAINTENANCE ACTION, START WORK 
ON THE NEXT MAINTENANCE ACTION IN THE Q, AND MOVE ONE THE REST 
UP ONE POSITION IN THE STACK 

r hr r r r ir r r r r r de de de ŷr ŷr ŷr ŷr de de de de de de de de r r r r r r ŷr ŷr ŵr ir ir kr vr vir ŷr ŷr ŷr ŷr ŵr ŷr ŵr ŵr fc ŵr ŵr ir ŷr ic ŷr r r r r r ŷr ŷr ŵr ŵr ŷr ŷr ŷr ŵr ŷr fc ŵr fc ŷr ŷr 1 1 1 š š š 3: ŷr ŷr ŷr ŷr r r or r 13 š 1 ŷr ŷr ŵr 
55, GO TO (26,56,46,56,66,/76,86,96),ELAWMT(L) 
36 CHKIN(2) = CHKINQ(1) 

TTMINQ(2) = TTMINQ(2) + (CLOCK - CHKINQ(L)) 

DO 37, I = 1,NAMTMT - 1 
CHKINO(I) z CHKINO(I*L) 
ELAWMT(I1) = ELAWMT(I+L) 

37 CONTINUE 


e de de de de de e de de de e e k k k kt k k k k k kt k k kk de de k k kt k k k k k k k k k K k k k k k k k Kk Kt k k k k k k k k k k k k k k k k kt k Ŵr k k k k k kt k k k k k k k k k k k k k kt k k k k K k ŵr ŷr jr ŵr ŵr ir K 


UPDATE LENGTH OF QUEUE 
defe dee defe fe dede dede fede dede e dee e dee kO K K ee e e e de k de D k Ñ e de e e e de Se e de Se e de e e e de de e de e de e de e dr dr dr dr e e e de e de t de de Se hr dr de e de de de Se Se de de ár dr dr dr dr dr dr dr dr de de de de e 
NAWTMT - NAWTMT - 1 
SVT - NORML(MSTOFO,SGMOF2) 
FEL(10) - CLOCK * SVT**(L/.33) 
RETURN 


*fe e fe de fe k k k k k k kt kt k k kt k k k k k ft k k k k k fr k kt k Kk k k k k k k k k k k k k k k k k k k k k k k k k k r k k k k k k k k k k k K k k k k k k k k k k k k k kk kt k k k k k k k k k k Kk k k k k K 


SCHEDULE NEXT BREAK-DOWN FOR 
Tr ik k ККК ККК De de Se e de De de Se de De De De de Se dr جل جلا جلا‎ de de e De e جل‎ de Se e e dee de de de de جلا‎ e e de e k K K k e de e dede dr dr dr dr e de de e e de de de dee d Se de de KO KOK KOK K 
38 IAT = EXPON(MITOF2) 

FEL(2) = CLOCK + IAT**([/.8&) 

RETURN 


"Òc fe e k r tr ir Yr ç k ŵr ile ŵr k k it ŵr k De k k k k k k k k K k k k k k k k k De k k k k k k k k k k k k k k K K k k k K k k k k k k De k k k k k k k k k k k k k K k k k k kk k k k k K K k K k K k k k k k K K K 


UPDATE CUMULATIVE STATISTICS 
3e de de de de RR Se Se АСЕ 
40 TREPTM(3) - TREPTM(3) * (CLOCK - CHKIN(3)) 

TNREP = TNREP + 1 

IF((CLOCK - CHKIN(3)).GE.3.0) F = F + 1 


"Èc `r fe e de de de e fe kt k k k k k kr de k De De k k k k k k k k k k k k k k k kk k k kt k k k k k k k k k k k k k k k k k k k k k k k k k kk k k k k k k k k k h K K k k K k k K K K k k k k K K k K K K K K K K K K 


CHECK CONDITION OF QUEUE 


Ve de "Òr de dede k k k k k k ft k k kt k kt ŵr Ŵr k kt k k K k k k k k k k k k k k k k k De k k k k Kk Ŵr k k k k k k k k k k k kk k k ŵr k k k k k k k t k k k k k k k k k k dr kt k k k k k K K k K K dr KOK K k K K K K K 


IF(NAWTMT).GE.1.0) GO TO 45 
FAI e de IOI e dk e e k k e e k Se de De e e de de De de de Se e e de Se k K e e e e e de de e de de e e ee de de de Se e جلا‎ de e de e e de Se de de de de de de de de de de جلا‎ de de de de de de Ho de ŵr W K K 
QUEUE EMPTY, SERVICE MAN BECOMES IDLE, NEXT DEPARTURE TIME 
EQUAL TO INFINITY 


49 


Wc cr ŵr îr îr cr ŵr k k r k cr ŵr de de e e Ve k K K K k k k k k k K k de cr ŵr ŵr r de k k K K ŵr ŵr rr ŵr ŵr ŷr ir ŵr ŷr K ŷr ŵr rr 0 K k kk ŵr ŵr 0 K îr k K ŵr ŵr ŷr k kk k k K de de K ŷr Ŵr ŵr ŷr ŵr rr ŵr ŵr de De âr ŵr k K K ŵr Yr œ 


NMTBUS - NMTBUS - 1 
FEL(11) = 1.0E+30 
GO TO 48 
Tc de dr Hc Üe de de Ve Ve de Ar Ae Ye КИЕ ŵr ŷr ŵr ŷr ŵr ŵr ŷr ŵr Ïe ŷr ŵr ŵr ŵr ŵr ŷr ŵr ŵ ŷr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr de ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr de Ез ХЕ ŵr âr ŵr ŵr ŷr ŵr ŷr ŵr de ŵr ŵr de ŵr ŵr ŵr r ŵr r ŵr ŷr Ŷr ŷr ŷr ŷr de ŷr de ŷr Ŷr ŷr Ŷr 
OUEUE CONTAINS AT LEAST ONE MORE MAINTENANCE ACTION, START WORK ON 
THE NEXT MAINTENANCE ACTION IN THE O, AND MOVE-THE REST UP ONE 
POSITION IN THE STACK 
ŵr Tc ŷr e ŷr ŵr КК КААК КЛИК ŷr ŵr ŵr ŵr ŷr ir ŵr ŷr ŷr ŵr ŵr ŵr ŷr ŵr ŷr ŵr r de ŵr rr ŵr ŷr ŵr ŷr ŷr ŵr r ŵr ŷr ŵr ŷr ŷr Ŵr ŵr ŵr ŷr r ŷr ŵr de ŵr ŵr ŷr de de ŷr ŷr ŷr ŵr ŵr ŷr ŵr ŵr ŷr de ŷr ŵr ŵr ŵr ŵr ŷr r ŷr We ŷr ŵr ŷr ŵr ŵr ŵr ŵr ŵr ŵr ŷr ŵr ŵr 
45 GO TO (26,36,46,56,66, 76, 86,96), ELAWMT(L) 
46 CHKINC3) - СНКІМО(1) 
TTMINO(3) = TTMINO(3) * (CLOCK - CHKINQ(l)) 
DO 47 1 s 1,NAWTMT - 1 
CHKINO(I) = CHKINQ(I*l) 
ELAWMT(I) 5 ELAMMT(I*L) 
47 CONTINUE 


Dc بار جار‎ Be de И ede fe eee eee ee eee ŷr hee ŷr de rc k K ŵr k ŵr fe ee ŵr K K Kt k r K ŷr ŵr k K ŵr k K kt k k K fe ee ee ee ŵr ŵr ŷr ŷr ŵr ŷr ŵr e eee ee efe r K k k K K k K r vr r r or r or o Ax 


UPDATE LENGTH OF QUEUE 
rdr de de He Ae de Ae de de de de de de de Ae de de de ŵr e fe e eee eee efe fe eee de de de de de de de de de de de de de De de de de de De de de de e he he he he he he he ke k e ee ee efe ŵr ŵr ŷr ŵr e efe ee ee ehe ee de ŵr ŷr Ŷr ŷr ŷr r 
NAWTMT = NAWTMT - 1 
SVT = NORML(MSTWS1,SGMWS1) 
FEL(11) = CLOCK + SVT**(L/.32) 
RETURN 


ŵr ŷr ŵr ŵr rir ŷr rir ir ŵr ŵr cr De ŵr ŵr ŵr ŷr ŷr ŵr ŵr ŵr ŵr ŷr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr r r De cc ŵr ŵr ŵr ŵr ŵr ŷr ŷr ŷr cc ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr de rr ŵr ŵr ŵr ŵr ŵr ŵr e e Ye yr r r cc rr rc ccc ŵr ŷr e de de de De De De ŷr ŷr ŷr ŷr ŷr ŵr ŵr de de Àe 


SCHEDULE NEXT BREAK-DOWN FOR WST1 


ŵr ŷr ŵr ŵr ŵr Ye De ir ir De ŷr ŵr r ŷr ŵr r de De ŵr ŵr rc ŵr ŵr ŷr ŵr ŵr ŷr rc ŷr ŷr ŷr e Âr ŷr ŵr ŵr ŷr ŵr ŵr ŵr rc ŵr ŷr Âr ŷr ŵr ŷr ŷr ŵr ŵr ŵr ŵr ŵr ŵr ŷr ŷr ŵr rr ŵr ŵr ŷr r ŵr ŵr ŷr rc ŵr ŷr De ŵr Âr ŷr ŵr ŵr ŷr rc ŵr ŷr ŵr de cr ŵr ŵr ŷr ŵr ŵr ŷr ŵr ŷr de ŵr ŵr de ŵr ŵr ŵr ŷr ŷr ŵr ŵr 


48 IAT - EXPON(MITWS1) 
FEL(3) = CLOCK + IAT**(l/.7) 
RETURN 


ŷr ŷr ŵr ŵr fe he ir e hne e ŷr De ŷr ŵr r ŵr ŵr ŵc ŵr ŷr ŵr ŵr ŷr ŵr ŵr ŷr ŵr ŵr ŷr ŷr ŵr r ŷr ŵr ŷr ŷr ŵr ŵr ŵr ŷr e ŵr ŷr ŷr ŵr De ŷr ŷr ŵr ŷr ŵr ŵr ŵr ŷr ŷr ŵr e ŷr ŵr de ŷr ŵr ŵr ŵr ŷr ŵr ŷr ŷr ŷr ŵr ŷr ŷr ŵr de de de de de de ŵr de de de de de ŵr rr ŵr rr ŵr rr ŷr rr ŵr ŵr e 


UPDATE CUMULATIVE STATISTICS 
x Ë We W dr dr Ach Sir do ы ы âr Wr Wr yr Wn rrr ЛҮ dh dh dh dh dh dh dr dh Be Be Br Sr Sr dèr Dr de De Be De ie Be dèr de Se Se cde de eye tolle de d de de de Wr k k k W 
50 TREPTM(4) = TREPTM(4) + (CLOCK - CHKIN(4)) 
TNREP = TNREP + 1 
IF((CLOCK - CHKIN(4)).GE.3.0) F = F + 1 


St جار جار‎ de ŷr e Wr De de e de de Ve Ve De diede Ve de de Ye De de ŵr ŷr ŵr ŷr ŵr ŵr ŵr ŷr ŵr ŵr e r ŵr Yr ŵr yr ŵr e de de ŵr cr ŵr ŵr ŷr de Ve e ŵr c ŵr Ye ŵr de de de de de de ŷr ŵr Âr e ŷr ŵr Yr ŵr Yr ŵr Yr r ŵr ŵr rr rc ŵr ŷr rr ŵr ŵr ŷr de ŵr جل‎ ŵr cr de de ŵr de Àe r o r 


CHECK CONDITION OF OUEUE 


He He c e ŷr ŷr cr ŷr ŵr ŵr ŷr ŵr ŵr ŷr ŵr ŵr ŵr ŷr ŵr ŷr ŷr ŵr ŷr ŷr ŵr rr ŵr ŵr ŷr r ŵr ŵr ŷr ŵr ŷr ŷr ŵr ŷr ŷr ŷr ŵr rr ŵr ŵr ŵr ŵr ŵr ŷr ŷr ŷr ŵr de de de de ŵr ŷr ŵr ŵr de ŵr De ŵr ŵr ŵr ŵr r ŵr ŵr ŵr ŵr ŵr ŵr ŷr ŵr ŵr ŵr ŵr ŷr ŷr ŵr ŷr ŵr ŵr ŷr ŵr ŵr e de Ye ŵr ŷr de ŵr ŵr ŵr ŵr Hœ 


IF(NAWTMT.GE.1.0) GO TO 55 
De de ide hn e et ee ee e ŷr ŵr de de ŷr ŷr ŷr ŵr ŵr ŷr ŷr ŵr ŷr ŷr ŵr ŷr ŷr ŵr ŷr ŵr ŷr ŵr ŷr ŷr r ŵr ŷr ŷr ŵr de e ŵr ŷr ŷr ŵr ŷr ŷr ŷr ŵr ŷr ŵr ŵr ŷr ŷr ŵr ŵr ŵr ŵr ŵr r ŵr ŵr ŵr ŵr ŷr ŷr r ŵr ŵr ŵr de ŵr ŵr ŵr ŵr ŷr ŷr ŷr ŵr ŵr de ŷr ŵr ŵr ŵr ŷr ŵr ŷr îr Me ŵr ŵr r ŵr x 
OUEUE EMPTY, SERVICE MAN BECOMES IDLE, NEXT DEPARTURE TIME 
EQUAL TO INFINITY 
Wc ir ŵr ŵr ŷr de Ve de de ŷr ŷr Ve ŷr ŷr ŵr ŷr Ve de de ŷr ŷr ŷr ŵr ŵr ir. ŵr ŵr ŷr ŵr ŷr rc ŷr ŵr ŷr ŵr ir ŷr ŷr rc ŵr rc ŷr Ŷr ir ŷr ŷr rc ŷr ŷr ŵr ŷr ŵr ŷr fe ŵr ŵr ŵr ŷr ŵr ŵr ŵr ŵr de r ŵr de ŵr ŷr ŵr ŵr r ŵr de de ŵr ŵr ŵr جل‎ ŵr de ŵr fe ŷr ŵr ŵr de de de ŵr de ŵr ŵr de de de de ŵr n Àe 
NMTBUS = NMTBUS - 1 
FEL(12) = 1.0E+30 
GO TO 58 
Yr ŵr ŵr ŷr K k kt k kk cr dx De de De k k îr ŷr ŷr K ŷr ŷr ŷr K K ŷr ŷr K ŵr ŷr ŷr ŵr ŷr ŷr ŷr ŷr ŷr ŵr ŷr K k K k ŵr ŷr K ŷr ŵr ŷr ŵr ŷr ŷr ŵr ŷr ŷr ŵr K K ŵr ŷr ŷr ŵr âr r ŷr k K t ŵr ŵr ŵr r ŵr ŵr ŵr ŷr ŷr ŵr ŵr r ŵr ŷr ŵr ŷr ŵr ŷr ŵr de ŵr r ŵr ŷr ŵr rc ŵr ŷr n x 
QUEUE CONTAINS AT LEAST ONE MORE MAINTENANCE ACTION, START WORK 
ON THE NEXT MAINTENANCE ACTION IN THE Q, AND MOVE ONE THE REST 
UP ONE POSITION IN THE STACK 
yr rr ŵr ŷr ŵr ŵrŵr ŵr ŷr ŷr ŷr ŷr ŵr r ŷr ŷr ŵr ŷr ŷr ŷr ŷr ŵr de ŵr ŷr ŷr r ŵr rc ŷr ŷr ŷr ŵr ŵr ŵr r ŷr ŵr ŷr ŷr ŵr ŵr ŵr ŵr îr r yr yr îr ŷr ŵr r ŷr ŷr ŵr e ŵr ŷr ŷr ŵr ŵr ŷr ŵr جل‎ ŵr ŵr r ŵr ŷr ŷr ŵr ŵr ŷr ŷr de ŵr ŵr ŵr ŵr KK K ŵr rr ŷr ŷr ŵr ŵr ŵr ŵr ŵr r ŵr ŵr ŵr ŵr ŵr * ŵ 
55 GO TO (26,36,46,56,66,76,86,96),ELAWMT(L) 
56 — CHKIN(4) = CHKINQ(1) 
TTMINQ(4) = TTMINQ(4) + (CLOCK - CHKINQ(L)) 
DO 57, I = 1,NAWTMT - 1 
CHKINQ(I) » CHKINQ(I*l) 
ELAWMT(I) - ELAWMT(I*L) 
57 CONTINUE 


T irit k fe t kit t جل‎ Dde ir k k k K k K ŵr rr de de de جل‎ K ŵr de e e K cr rr k k kr e e e ee جل‎ k cr k k k k k k k r k k e rr rr rr k K k K k cr de k k cr kt e K k k k ŵr ŷr ŵr ŷr ŵr ŷr ŵr k ŵr ŵr ŷr ŵr K ŵr ŷr ŵr ŵr K k K k e 


UPDATE LENGTH OF QUEUE 
LLATAI 

NAWTMT = NAWTMT - 1 

SVT = NORML(MSTWS2, SGMWS2) 

FEL(12) = CLOCK + SVT**(L/.32) 

RETURN 


He Tc èc e e e d e e ee ee e e k r t K K K K k k k k k r k k e ŷr ŷr de k K K k k k r de de de de de r cr k k k K k K k k k k k k K k de k K r k k k Ŵr K K K k K k k k ŷr k K K K ŵr ŵr r ŵr k K k k k k ŵr ŵr ŷr ŷr ŵr 


SCHEDULE NEXT BREAK-DOWN FOR WST2 


3t Se dc de Ye Te L de De Ue k k k k k Ŵr Ŵr K k k kt k t k Ŵr k k k K K k k k k t de ŵr k de K k K k K K k h k K e k k t k t k جل‎ k t k k k k k k k k k k k kk ŵr ŷr k K kt K k k k k k cir Wr k K k k k k k k k x 


50 


58 IAT = EXPON(MITWS2) 
FEL(4) = CLOCK * IAT**(l/.7) 
RETURN 


d de fe e ee dede ee de cf ee dece dne eee rite eee eee fee eee eode eee ee eee e ee eee o eee e e eee ee ee e e f de k KK K k KE œ 


UPDATE CUMULATIVE STATISTICS 
EREREAAREE e de de Se de de de de de de de de de dede eee dede dede de dede dede dede de de de de V de d de K dede de de 3 K d d R K K d d K d K K K d d K d d K í K K K d K K Re e OR GO 
60 TREPTM(5) = TREPTM(5) + (CLOCK - CHKIN(5)) 

TNREP - TNREP * 1 

VECCCLOCK - CHEIN(DOO. GE. 5.0) F = F + 1 


как} ССК КККК ККИ КККК КК ИККО X d d H КТЖ ККЕ КЕК КЕ ИКЕ ЖЕКЕ * * * * * * * * * * * * 


CHECK CONDITION OF QUEUE 


ЖЕКИ ОИЕ EE X À À $c * * H H Hœ HŒ H H 39: H H H Wœ d КККК ЇКЇ КККК КККК КККК КККК КККК ККК КККК ККК КККК КККК КККК КККК ККК КК ККК КК 


IF(NAWTMT).GE.1.0) GO TO 65 
kA he e fee fedet dede dede eee dehet eee c ee ee fe eee fe eee e eee eee ee eee ere e eee eee ee eee eee fe d d 
QUEUE EMPTY, SERVICE MAN BECOMES IDLE, NEXT DEPARTURE TIME 
EQUAL TO INFINITY 
Tr Tx Ie A 337 e de Ve die Ve Ve De de de Òe e LA AAA AAA AAA AA AA AAA AAA AAA AAA AAA AAA AAA AAA AAA AAA AA AAA AAA AAA AAA AA AAA AA AAA AAA AAA AAA AAA 
NMTBUS = NMTBUS - 1 
FEL(13) = 1.0E+30 
GO TO 68 
* k k kk k k fr fr fr k Me fr k k k k k k k k k k k N k k k k k k k t k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k kk k k kk k kk k k k k k k k k t k kk k k k 
QUEUE CONTAINS AT LEAST ONE MORE MAINTENANCE ACTION, START WORK ON 
THE NEXT MAINTENANCE ACTION IN THE O, AND MOVE-THE REST UP ONE 
POSITION IN THE STACK 
fe e k k k k k k k ft fr kr k k k k kk k k k k k k k k ff k k k f ff k k Me k Me k k ŵr k k Me ŷr ff k t k ŵr k ôr ŵr ŷr k ŷr ŷr ŷr k k kt ŷr ŷr de ŷr ŵr ŷr ôr ŷr ŷr k k k fk k k ŷr k k t f k f ff k ff ft k k f f k k k k k k k k k k k 
65 GO TO (26,36,46,56,66, 76,86,96) ,ELAWMT( |) 
66  CHKIN(S) = CHKINQ(1) 
TTMINO(S) = TTMINQ(S) + (CLOCK - CHKINO(L)) 
DO 67 1 = 1,NAWTMT - 1 


CHKINQ(I) = CHKINQCI+1) 
ELAWMT(I) - ELAWWT(I*L) 
67 CONTINUE 


* kk k k k k k fr k k k k k k k k k ft k kk k ft k k k k k k k k k k t k k k k k k k k k k k k k k k k kk l k k k k k k k ft k k k k k k k k k k kt k k k k k k k k k kk k k k k k kk k k k k k k k k k k ir 


UPDATE LENGTH OF QUEUE 
e de k k k k k k k fr k fr k k k k k k k k k k k k k k k k k k k k k k k k k k k k kk k k k k k k k k k k k k k k k k k k k k k kk k k kk kk kk kk 39: kk kk kk k k k k k kk k kk k k k kk k kk k 
NAWTMT = NAWTMT - 1 
SVT = NORML(MSTPT1,SGMPT1) 
FEL(13) = CLOCK + SVT**(l/.5) 
RETURN 


T fr k k k k fk k k k k k k k k k t k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k ff k k k k k t k k k k k 


SCHEDULE NEXT BREAK-DOWN FOR PTT1 
Pe e Be He Hh De D K Me de Me de de de de De Se M de de M de de de de De de de de de de Е 
68 IAT = EXPON(MITPT1) 

FEL(S) = CLOCK + IAT**(l/.77) 

RETURN 


fefe defe he de fe k k k k Tr k fr fr k t k k k k k K k k k k k k k k k k k k k k k k k k k k k k k k k k k k f k k k f k k f k k k t ff k k t k k k k k k k t k k t ft k k k k k t ft f f ft ft k k k k k k k k k ft k 


UPDATE CUMULATIVE STATISTICS 
KHKKKKKKKKHHKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKHKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKKK 
70 TREPTM(6) = TREPTM(6) + (CLOCK - CHKIN(6)) 
TNREP = TNREP + 1 
IF((CLOCK - CHKIN(6)).GE.3.0) F = F + 1 


+ t de k kt k k k k k ft t k fr t k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k kt f k f k t k k t f k ŵr ŵr k k k k k k k k k k k k f k ŷr ŷr ŵr ŷr ŵr ŷr k k k k ft k k k k k k œ 


CHECK CONDITION OF QUEUE 


k ftrt fr k k k k k fr + Me k k k k k k k k k k k k k k k kk k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k k 


IF(NAWTMT.GE.1.0) GO TO 75 
ŵr fefe de de KE dee N H De HK HK TH OK fte oed eee eee hr H d d d d fee eoe toon feo ette oe ee oe eee oed eoe n d d A 
QUEUE EMPTY, SERVICE MAN BECOMES IDLE, NEXT DEPARTURE TIME 
EQUAL TO INFINITY 
جار جا‎ fee de de efe de de dee dede dee dede k k k k ft De Me ŷr ŷr k h k t k Me ft De eite tee eee de Me Me Me Me Me de Me de de de de de de Me de de Me Me de de de Me de de de k k k ŵr k ŵr ŵr ic k k k ŵr ŵr ff k k k k k ft k k k k k k k K 
FEL(14) = 1.0E+30 
GO TO 78 
keKktkKkkekkk&kkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkkk&kkkkkk&kkkkkhkkkkkkkkkk&kkkkkkk&k&kkkkkkkkkkkkkikk 
QUEUE CONTAINS AT LEAST ONE MORE MAINTENANCE ACTION, START WORK 
ON THE NEXT MAINTENANCE ACTION IN THE Q, AND MOVE ONE THE REST 
UP ONE POSITION IN THE STACK 
* de fe k k k kk k k k kt k k k k kk kk kk k k k k k k k k k k k k k k k k k k k k k k k k k k k ff k k k k k t k f f k k k k f k k k k k k k k k k k kk ft ff k k k k k k k k k f k k k k k k t k k % K 
75 GO TO (26,36,46,56, 66, 76, 86,96), ELAWMT(L) 
76 CHKIN(6) * CHKINQ(1) 
TTMINQ(6) = TTMINQ(6) + (CLOCK - CHKINQ(l)) 
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DO 77, 1 = 1,NAWTNT - 1 
CHKINO(I) = CHKINQ(I+L) 
ELAWMT(1) = ELAWMT(I+L) 

77 CONTINUE 


fikk kiki i Yr KE e kr K k k X k k kk K k k k k k k K k Kk k k k k K k k k k kk k k k K K k k k k k K K K k k k k k k k k K k k k K k K K k K k k k k k K k R K k R K K R K k R k k k k K 


UPDATE LENGTH OF QUEUE 
x à i re vc D dc Tr hr Ar АСААВ e e e 
NAWTMT = NAWTMT - 1 
SVT = NORML(MSTPT2 , SGMPT2) 
FEL(14) z CLOCK * SVT**(l/.5) 
RETURN 


fk k k k'iri k Yr ККА КИК k k K X 0 k KK K k k K K K K K Kk K K k kk K k k k K K K K k k K K k k K K K K K K k k kk k k k k k K k k K K K K K ŵr k k K K K K ŵr ŵr KOR RO K KOK R K K R K K K K Ж 


SCHEDULE NEXT BREAK-DOWN FOR PTT2 
«KAAKA AAA AAA Ñ K KA AKKA AKKA de e IR dh de De e dh De de AKKA KAKAK KA KAKAK KAAKA KAKA KAK KA AKA KAK KAA AKKA AKKA AK KAA AAA 
78 IAT z EXPON(MITPT2) 
FEL(6) = CLOCK * IAT**(l/.77) 
RETURN 


Sede fee de dede de ee ee ee ee eee ee eee جل :ل جل‎ dee eee eee e eee eee ee ee eee ce eee eee e eee ee eee ŵr ŵr ŵr rrr rr ŵr ŵr rr ŵr ŵr fe ŵr dee 


UPDATE CUMULATIVE STATISTICS 
Www rd rw dr rr rir rr rrr rrr RR RR RR RR RR RR RH RH eee ee eee oe de e de e e e 
80 TREPTM(7) = TREPTM(7) + (CLOCK - CHKIN(7)) 

TNREP = TNREP + 1 

IF((CLOCK - CHKIN(7)).GE.3.0)F - F * 1 


t de dede dede de deed dece e defe * * * * * * * * E * * * îr * * ** * ** * * * * * * * * * ŵ * * * * * * * * * ** * * * * * * *c* ККЕ ККЕ КЕ جل‎ r 97 r r جل جل جل أ جل‎ r ŵr * de dede dede * * * 


CHECK CONDITION OF OUEUE 


E E ir r r rk Î r Û kk Û Û r O Û Û Û i r r E r E r o r r r Û E Û r r r r kr r kr r E r r r r r r r r r r r r r r kk r r r r r Fk ok Fk kk 


IF(NAWTMT).GE.1.0) GO TO 85 
Yk Tk He He de he he de de He ir de HU HX X ŵr de de he de  * * * He ŵ de * * * * * * * * * ŵr ŵr he he de ŵ * ŵ * ŵ * * * ŵr * ŵ * * * * * * de جل‎ * * * ŵ * * ŵ de ŵ * ŵ de جل‎ ŵ * de ŵ * de de * ŵ * * * ŵ * * * * * ŵ * * * ŵ * * 1 
OUEUE EMPTY, SERVICE MAN BECOMES IDLE, NEXT DEPARTURE TIME 
EQUAL TO INFINITY 
XXX ÈX X de dedere de eee e eee ee eee ŵ î r ŵ * * * * $ ŵ ŵ r ŵ ŵ ŵr ŵ ŵ * * ŵ * ŵ ŵ * * ŵ * * *  * * ج11‎ * ŵ جل * جل‎ ŵ ŵ ŵ ŵ ŵ * * ŵ ŷ جه جل جل جل‎ ŵ ŵ ŵ ŵ جل جل جل جل جل * جل جل‎ ŵ ŷ UUU n d 
NMTBUS = NMTBUS - 1 
FEL(15) = 1.0E+30 
GO TO 88 
dfe dede de de de dede جلو جاو‎ ŵr ic ŵr ŵr ŵr ŵr ŵr جا‎ ŵr * * ŵr * ŵr * ŵr ŵr ŵr ŵr ŵr ŵr ŵr îr îr ŵr îr ŵr ŵr îr îr ŵr ŵr îr ŵr ŵr ŵr ŵr îr ŵr îr îr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr îr ŵr * * r * ŵr * ŵr * ŵr ŵr ŵ ج11‎ * ŵr * îr * * * * * * ŵ ŵ * * * * * * 
OUEUE CONTAINS AT LEAST ONE MORE MAINTENANCE ACTION, START WORK ON 
THE NEXT MAINTENANCE ACTION IN THE O, AND MOVE-THE REST UP ONE 
POSITION IN THE STACK 
ŵ rc ŵ КККК КК Кї ККК îr ŵ * ŵ ŵ îr ŷr ŵr îr îr ŵr ŵr ŵr ŵr ŵ îr ŵr ŵr îr ŵr ŵ ŷr ŵ ŵr ŵr îr ŵr ŵ îr ŵ îr ŵ îr îr ŵ ŵr * ŵr * * * ŵr ŵ * * ŵr ŵr * * * * * RK îr ŵ ŵ * ŵ * îr ŵ îr îr ŵ ŵr * ŵ îr ŵr ŵr ŵr ŵr d * n x 
85 GO TO (26,36,46,56,66,76,86,96),ELAWMTCL) 
86  CHKIN(7) = CHKINQ(1) 
TTMINQ(7) = TTMINQ(7) + (CLOCK - CHKINQ(L)) 
DO 87 1 = 1,NAWTMT - 1 
CHKINQ(I) = CHKINQ(I*l) 
ELAMMT(I) s ELAWMT(I*l) 
87 CONTINUE 


tfe fe efe fece eee ce dede جل جل‎ fe dee e ec îr ŵr îr ŵr ŵr ce ce eee ee eee ee ode ee eee e e nee ee eee eee eec cce eee cec ce e cce e ee c efe fefe * * * * * * 


UPDATE LENGTH OF QUEUE 
krik er der ŵr droe dr or der dh dr dh dh dr dr dr fr tr tr drdd ddr dor dr dr dr dr ŵr dr dr dr dr dr dr dd dr Be Pe Se Se Te dr de e Te Te Te Te Te dr dr dr dr dh dh dir e He Se de De d de de de e de de dr de de de de e e e e 
NAWTMT = NAWTMT - 1 
SVT = NORML(MSTPC1,SGMPC1) 
FEL(15) = CLOCK + SVT**(L/. 14) 
RETURN 


fikk k'ik kt k k k'ik; K k k kk k k K k K K K K KK kk kk k K k KK K kk kk k k k K k K K k k k k k k k k 0 k k k K k k K K k kk kk k k k k K k k k K K K K k k k k k K KK KK k k K K k KKK KKK KKK 


SCHEDULE NEXT BREAK-DOWN FOR PTT1 
КАЖАА ЕАТС ЕАК НАВОИ НИКИТА dee d io 
88 IAT = EXPON(MITPC1) 

FEL(7) = CLOCK + 1AT**(1/.73) 

RETURN 


dese fee dede fede dee dede dede fede de dee e k k k K kk k kk k k k k K k k k k k k k k k k k k kk k k k k k k k kk k k k k É Kk k k k k k k k k k k kk k K k k k k kk k k k kk k k k k k k k k k k k k k k k k k K: 


UPDATE CUMULATIVE STATISTICS 
кезе а ка КИЕК ЕК К ЕЕ КК КАЕ К Ж ht dr RR rr dr dr dr dr dr dr dr dr dr dr der e e hr dh dh Ár dr dh جلا جلا‎ dr dh dr dh de dr dh dh dhd dh dr d dr dh De de de جلا‎ Èc de He de OR ŵr ŵr e 
90 TREPTM(8) = TREPTM(8) * (CLOCK - CHKIN(8)) 

TNREP = TNREP + 1 

IF((CLOCK - CHKIN(8)).GE.3.0) F = F + 1 


TKK KKR AAR KK KK LK جف جل جل :ل جف جل جل 117 ج11 جل جف‎ LLL LLL LKK KL جف جف‎ he eee che de fe ee cde e eee eee r r r dede o dede de hee de efe de r r r r r r r r eh 


CHECK CONDITION OF QUEUE 


c Tc Ic dc Ie k kk k'ir kk k k k k r k f e ccc ecce ee k k k k k k k k k k Kk kk kk kk k k kk k k k k k k k k k kk k kk kk k k kk k k k k k k e eec ch ne ede e e k k K k k K k Wœ Wœ e e d d Wœ 


IF(NAWTMT.GE.1.0) GO TO 95 
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Wr ir ŷr Vc de ŵr vk vir r e ee ee e e ŵr ŵr K ŵr ŵr ŵr ŵr ŵr ŵr Ve De de de De Òe de K K k k k k k k k k k k k k K K K K k k K K k k k de k K K KK Kk k Kk de K k k k k k k k k k k k k K K k k De k k k K k K k k k k Kk K K K 


QUEUE EMPTY, SERVICE MAN BECOMES IDLE, NEXT DEPARTURE TIME 
EQUAL TO INFINITY 
«AA AAA KAA AAA AAA А e e De Tie Tie De de De de e Se e De e de de de e De e De De de de de de de de de De de Ik 
FEL(16) = 1.0E+30 
GO TO 98 
«кукук ККК КК ИККЕ III IOI De علا‎ 4+ de d de d جل جل جف جف‎ De de de 9> جا جا جا جا جا‎ De de De d de Se Se جف‎ De de جلا جلا جلا جلا جلا‎ de de de جلا جلا جلا‎ de de de de de de tt tick 
QUEUE CONTAINS AT LEAST ONE MORE MAINTENANCE ACTION, START WORK 
ON THE NEXT MAINTENANCE ACTION IN THE Q, AND MOVE ONE THE REST 
UP ONE POSITION IN THE STACK 
dA AKA AAA AA AA A A AAA K KK AAA КККК ИК ИК ККК ККК КОК ОИ КККК КОК КК ККЕ КЕ De de de de de de De de De de e de de de de de de de ee ode de de CE e d d R R R e de 
95 GO TO (26,36,46,56,66,76,86,96),ELAWMT( L ) 
96 СНКІМ(8) г СНКІМО(1) 
TTMINO(8) - TTMINO(8) + (CLOCK - CHKINQ(l)) 
DO 97, I = 1,NAWTMT - 1 
CHKINO(I) -z CHKINOC(I*L) 
ELAWMT(I) = ELAWMT(I*L) 
97 CONTINUE 


Sce fede ôr ŷr ŷr ŵr ŵr ŷr ŷr ŷr ŷr ŷr ŷr ŵr ŵr ŵc ŷc ŷr ŷr ŷr c ŵr ŵr ŵr ŵr ŵ * ŵr ŷr ŷr ŵr ŵ ŷr De ŷr ŵr ŵr ŵr ŷr r r ŵr de de ŵ ŵr Ŷr de de Ŵr ŵr ŵr ŵr ŷr ŵr ŵr Ŵr Ve ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŷr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŷr ŵr ŵr kkk k 


UPDATE LENGTH OF OUEUE 
rc de d e e e e e e e e e e e e e de c d e e e e ee eee dede ee He e e e IIIS IIIS ISIS OI KA KA A 
NAWTMT = NAWTMT - 1 
SVT = NORML(MSTPC2,SGMPC2) 
FEL(16) = CLOCK + SVT**(1/.14) 
RETURN 


Yr ŷr ŷr ŷr ŷr ŷr ŵr ŵr rr ŷr ŷr ŵr ŷr ŵr ŷr ŵr ŵr ŵr ŵr ŵr ŵr ŵr Tr fr Tr Tr Tr r r r K ŵr ŷr ŵr ŷr Yr ŵr A ŷr ŵr ŷr ŷr ŵ ŵr k k ŷr ŵ ŷr ŵr ŷr ŷr ŵr ŵr ŷr ŵ ŵ ŵr ŵr ŵ ŷr ŵ ŵr ŷr ŷr ŵ ŷr ŵr ŵ ŷr ŵ ŵ ŷr ŵ ŵ ŷr ŵ ŵr ŷr ŵ ŷr ŵr ŵr ŷr ŵr ŷr ŵr de ŷr ŵr ŵr ŵr ŵr ŵ A ŵr ŷr ŵr 


SCHEDULE NEXT BREAK-DOWN FOR PCT2 
KAAKA KKK rr rrr rr rr rrr ŵr ŵr rr rr rr rr dr ŵr ŵr ŵr rr rrr rr rr dor de De de e de De de de de ddr dr dr dr dr dr dr rd de جل‎ e He dede de dede ke ke de de kekeke deke ke dek dek kekeke kekeke لاجلا جلا‎ 
98 IAT = EXPON(MITPC2) 

FEL(8) = CLOCK * IAT**(l/.73) 

RETURN 

END 


ŷc ŷr ŷr ŷr ŵ ŵ ŷr ŷr ŵ ŷr ŷr ŷr ŷr ŵr ŵr ŵr ŵr ŷr c ŷr ŵr ŵr ŵr ŵr ŵr ŷr ŵr ŷr ŵr ŵr ŵr ЖК De De De De de e de de de e de de КК КККК ККК КК de Òe de de ŷr Ve جلا جل‎ Ve de de جل‎ de Ve e de de de de de yr ye vk ŵr de ŵr ŷr ŵr e Ŵr ŵr Ŵr ŵr ŵ ŵr ŵr ŵr ŵ ŷr ŵ ŷr ŵr Hœ 


SUBROUTINE RPTGEN 


c c c К de De de e De De de Òe ŷr ŵr De De ŵr ŵr ŵr ŷr ŵ ŷr ŵr ŷr ŵ ŷr ŵr ŷr ŵr ŵr ŵr de ŵr ŵr Wc ŵr ŷr جل‎ wr de d de De de e ŵr ŵr Wr Wr Wc de d Ve * К de de de de de de e de de de de de de de d de d ŷr de Ye ŵr جل جل‎ ŵr ŵr ŵr ŵ ŷr r vk yr ŵr * de œ 


GENERATE REPORT OF SIMULATION RESULTS 

ŷc ŷr c ŵ ŷr ŵr ŷr ŵr ŷr ŷr ŷr ŷr ŵr ŵ ŷr ŷr ŷr ŷr ŵr ŵr ŷr ŷr ŵr ŵr ŵr ŵr ŵ ŷr ŵ ŵr ŵr ŷr ŷr ŵr ŵr Ек ŵ ŷr ŷr ŵ ŷr ŵr ŵ ŷr ŵ ŵr ŵr ŵ ŵr * * ŵr ŵ ŷr ŵ ŷr ŵr ŵ de ŵ r ŵ ŵr ŵ ŷr ŵr Ей de ŵ ŵr * ŷr ŵ ŵ ŷr Ŵ ŵ ŷr Ŷ КЕ ŷr Œ * ŵ ŵ 
INTEGER F,ELAWMT(100),TNREP 
REAL MITOF1,MITOF2,MITWS1,MITWS2,MITPT1,MITPT2,MITPC1,MITPC2, 
1 MSTOF1,MSTOF2,MSTWS1,MSTWS2,MSTPT1,MSTPT2,MSTPC1,MSTPC2, 
2 FEL(17),CHKIN(8), CHKINQ( 100), NORML, TREPTM(8) , TTMINQ(8) 
COMMON /SIM/ MITOF1,MITOF2,MITWS1,MITWS2,MITPT1,MITPT2,MITPCL, 
1 MITPC2,MSTOF1,MSTOF2,MSTWSL,MSTWS2,MSTPT1,MSTPT2,MSTPC1,MSTPC2, 
2 SGMOF1,SGMOF2, SGMWS1,SGMWS2, SGMPT1,SGMPT2, SGMPC1, SGMPC2, NAWTMT, 
3 NMTBUS,CHKIN,CHKINQ,TREPTM, MQ, F, TNREP, NMTPER,ELAWMT, TTMINQ 
COMMON /TIMEKP/ CLOCK, IMEVT,NUMEVS, FEL 


ŵ ŷr ŵr ŵr ŵr ŷr Wr ŵ ŵ ŷr ŵr ŵr ŵr ŵ ŷr ŷr ŵr ŵr ŷr ŷr ŷr K K K k k ŵr de ŵr ŷr ŵr ŵr ŵr ŵr ŵ ŵr ŵr ŵr ŷr ŵr ŵr ŷr ŵr k k ŵr ŵr e fe ŷr ŵr ŵr ŷr ŵr ŵ ŵr ŷr ŵ ŷr ŵr ŵr ŷr ŵr ŵr ŷr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵ ŷr ŵ ŵ ŷr ŵ k Ŵr ŵr ŵr ŵ ŷr ŵr ŵr ŷr ŵr ŵ ŷr ŵr ŵ í ŵ * * 


COMPUTE SUMMARY STATISTICS 


xe he e e fe De Ye De De Ve De De De de He Ve de de Ve de de vir de de De Ve de KL KEL LE e de de Ye r Ye Ve Ye e de r Ye e Ye e Ve Ye Ve de De Ve fe K KEL EKLE LELE KL KELL KL bK KKL bK bL KEL hL LL bb KA 


ОРКОҮ1 - (1-((TTMINO(1) * TREPTM(1)) / CLOCK)) * 100 
OPRDY2 - (1-((TTMINO(2) * TREPTM(2)) / CLOCK)) * 100 
OPRDY3 - (1-((TTMINO(3) + TREPTM(3)) / CLOCK)) * 100 
OPRDY4 - (1-((TTMINO(&) + TREPTM(4)) / CLOCK)) * 100 
ОРКОҮ5 - (1-((TTMINO(5) + TREPTM(5)) / CLOCK)) * 100 
OPRDY6 - (1-((TTMINO(6) + TREPTM(6)) / CLOCK)) * 100 
OPRDY7 - (1-((TTMINO(7) * TREPTM(7)) / CLOCK)) * 100 


OPRDY8 - (1-((TTMINO(8) * TREPTM(8)) / CLOCK)) * 100 


fk k K K Tr k rk fr r k Ŵr ŷr rr ŷr ŷr ŷr ŵr K ŵr k k k K ŵr Yr k e ŵr Yr De De de De k k e r k k k t Ye De t ie K Ж r ŷr ŵr ŵr ŵr ŵr ŵr ŵr ŷr ŵr ŵr ŷr r ke K fr fr ŵr ŵr rr ŵr ir ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr de k Ye x 


WRITE OUTPUT 


fkr k k rc ŵr rir t ŵr lr ŵr Yr k k fr جل جل‎ ŷr k KK جل‎ ŵr جل‎ ŵr K ŵr ŵr ŷr Yr t ŷr جل‎ ŵr جل جل‎ ŵr ŵr جل‎ K جل جل‎ rr ir rr rr ŷr ŷr ŷr ŵr ŵr ŷr ŵr ŵr e 0 k k k k k k kk K Wr K K O k k K K K K K k K k R OR R KOR R R KOR R K R RO R KO K Xx 


WRITE(6,10) 
10 FORMAT(6X, 'SIMULATION MODEL OF AIRCRAFT SIMULATOR MAINTENANCE', 

1///) 

WRITE(6,11) MITOF1,MITOF2,MITWS1,MITWS2,MITPT1,MITPT2,MITPCI, 

\MITPC2,MSTOFL ,MSTOF2,MSTWS1,MSTWS2,MSTPT1,MSTPT2,MSTPCL ,MSTPC2 
11 FORMAT(1X, ‘LAMBDA USED FOR OFT1 FAILURE TIMES ' F10.2 

1 /1X,'LAMBDA USED FOR OFT2 FAILURE TIMES 1, F10.2 

2 /1X,'LAMBDA USED FOR WSTI FAILURE TIMES ',F10.2 

3 /1X,'LAMBDA USED FOR WST2 FAILURE TIMES ' F10.2 
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4 /1X, "LAMBDA USED FOR PTT1 FAILURE TIMES ل‎ 
5 /1X,'LAMBDA USED FOR PTT2 FAILURE TIMES ',Е10.2 
6 /1X,'LAMBDA USED FOR PTCI FAILURE TIMES ',F10.2 
7 /1X,'LAMBDA USED FOR PTC2 FAILURE TIMES ', F10.2 
8 /1X, "MEAN SERVICE TIME FOR OFT1 REPAIR ', FlO.2 
9 /1X,'MEAN SERVICE TIME FOR OFT2 REPAIR ',FLO.2 
A /1X,'MEAN SERVICE TIME FOR WSTI REPAIR ',FlO.2 
B /1X,'MEAN SERVICE TIME FOR WST2 REPAIR ',FlO.2 
С /1X,'MEAN SERVICE TIME FOR PTTI REPAIR ' F10.2 
D /1X,'MEAN SERVICE TIME FOR PTT2 REPAIR ' F10.2 
E /1X, "MEAN SERVICE TIME FOR PTC1 REPAIR ' , FL0.2 
F /1X,'MEAN SERVICE TIME FOR PTC2 REPAIR ',F10.2) 
WRITE(6,12) SGMOF1,SGMOF2,SGMWS1,SGMWS2,SGMPTI,SGMPT2,SGMPC1, 
1SGMPC2 ,NMTPER 
12  FORMAT(1X, 'STANDARD DEVIATION FOR OFT1 REPAIR ',F10.2 
1 /1X, "STANDARD DEVIATION FOR OFT2 REPAIR ',F10.2 
2 /1X, 'STANDARD DEVIATION FOR WST1 REPAIR ',F10.2 
3 /1X,'STANDARD DEVIATION FOR WST2 REPAIR 1 F10.2 
4 /1X, "STANDARD DEVIATION FOR PTT1 REPAIR ', FlO.2 
5 /1X, "STANDARD DEVIATION FOR PTT2 REPAIR ',FLO.2 
6 /1X,'STANDARD DEVIATION FOR PTC1 REPAIR ',F10.2 
7 /1X, "STANDARD DEVIATION FOR PTC2 REPAIR ',FLO.2 
8 /1X, "NUMBER OF MAINTENANCE PERSONNEL ' 10.0) 
WRITE(6, 13)MQ, OPRDYI,OPRDY2 , OPRDYS, OPRDY4, OPRDY5 , OPRDYÓ, 
1OPRDY7, OPRDY8, F, CLOCK, TNREP 
13 FORMAT(LX, ‘MAXIMUM QUEUE LENGTH ' 110.0 

2 /1X,'OFTl OPERATIONAL READY RATE 1, F11.0, !X! 
3 /1X,'OFT2 OPERATIONAL READY RATE 1,F11.0,'X" 
4  /1X,'WST1 OPERATIONAL READY RATE ! a F11.0, ! X! 
5. /1X,'WST2 OPERATIONAL READY RATE ! F11.0, ! X! 
6 /1X,'PTT1 OPERATIONAL READY RATE 1 FLL.O, '%! 
7 /1X,'PTT2 OPERATIONAL READY RATE tr, FLL.O, 'X! 
8 /1X,'PTC1 OPERATIONAL READY RATE 1 FLL.O, ' X! 
9 /1X,'PTC2 OPERATIONAL READY RATE ELLO, 12! 
A /1X,'NUMBER OF MAINTENANCE ACTIONS THAT TOOK 1. 

B /1X,'MORE THAN THREE (3) HOURS TO CORRECT ',F10.2 

C /1X,"SIMULATION RUN TIME 1 F10.2 
D /1X,"NUMBER OF MAINTENANCE ACTIONS DURING SIMULATION '110.0) 
RETURN 

END 


Wr ŵr ŵr ŵr Ŵr Ŵr Ŵr yr جل‎ Ŵr Ŵr iir sr Ŵr sêr rr Ŵr r 5 5 3 5 5 5 5 5 3 5 D 3 5 5 3 5 3 5 5 3 5 3 3 5 d 5 5 5 5 5 h 5 5 3 5 h 5 5 3 5 5 5 D 8 5 8 5 5 5 5 d 8 5 5 D 5 5 8 5 5 5 5 5 3 3 8 5 5 5 8 5 8 5 3 5 5 5 8 3 5 2 d 
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Yr de dede de 2 2 2 2 4 $ x 2 rk x 4 2 X $. 4 2.2 2 ХКК ЖК r i ir r vir ir oir oir oir oir oir E oir oir oC oir oir oir oir oir oir dr cr ŵr ŵr ŵr e kr i oir i e e ice ee ŵr rr rr rr ŵr dr rr ŵr rr rr ŵr rrr rr ŵr rr ŵr rr ŵr Ж} 


EXPONENTIAL RANDOM VARIATE GENERATOR 
DEFINITION OF VARIABLES 


ALOG - Natural log function 

DSEED - Seed for random U(0,1) number generator 
EXPON - Exponential random variate 

FMEAN - Mean parameter of exponential distribution 
GGUBFS U(0,1) random number generator 

R - Random number 


rr de de e eddie e de en ee dede dede e he e eroe ee e e e ee e جا‎ ŵr e e e ee ee eee e eie ee її e ir ir Yr i ir Yr ir e ir ir ir i Yr ir ŵr rr ŵr vir oir oir ir Yr ir ir ir ir e e e le ir ir ir e ŵr ŵr de de e ŵr d d d ŵr 


FUNCTION EXPON (FMEAN) 
DOUBLE PRECISION DSEED 
DATA DSEED /43/ 
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GENERATE A U(0,1) RANDOM NUMBER 
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R=GGUBFS(DSEED ) 
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GENERATE AN EXPONENTIAL RANDOM VARIABLE WITH MEAN FMEAN 
3e de de e e eee ee de de e de e ee e ee e de de d e de dee De de de e dede de dede de e dede de ded de e dede dede deed de dede dede de e dee dee de dede dee dede dede doe dede dede de de de de de de de de de Ŵr OK 
EXPON--FMEAN*ALOG(R) 
RETURN 
END 
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Ŵr fc جاو‎ Ye ŵr hr ir ŵr ŵr ККК КККК К КК cr Ŵr ŵr cr dc dr ŵr câr dr ŵr ir dcr cr ir ôr ŵr cr iv cr de dece de dede de dr ŵr ôr dr ŵr ŵr dr ŵr ŵr dc Ŵr dr rr cr Ŵr e ir ŵr dr rr Ŵr cr rr ŵr rr ŵr rr ŵr ŵr dir ŵr ir dr ŵr rc ŵr ŷr ŵr ŵr ŵr ŵr ŵr rc ŵr ŵr 


NORMAL RANDOM VARIATE GENERATOR 
DEFINITION OF VARIABLES 


MEAN - Mean of normal distribution 

SIGMA - Standard deviation of normal distribution 
NSEED - Seed for U(0,1) random number generator 
NORML - Normal random variate 


K - Index to determine which method of computation to use 
PI - Pi 

RONE - A U(0,1) random number 

RTWO - Another U(0,1) random number 

ZONE - A N(0,1) random variable 


2TWO - Another N(0,1) random variable 
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FUNCTION NORML (MEAN,SIGMA) 

REAL MEAN, SIGMA, NORML 

DOUBLE PRECISION NSEED 

DATA K /0/, PI /3,14159/, NSEED /58/ 


fe de de e de e de e de de de dee ede de ede de de kt k k k k k k k kk k kk kk Kk k k k k k kk kk k k k k k k k k k kk kk k k k k e k de kk Kk k k k k k k k k k k k k k k k k kk k k k k k K kk k k kk k k k k k k k k k k ŵr ŵr ŵr K: 


CHECK TO SEE WHICH N(O,1) RANDOM VARIABLE TO USE 


Ve de de de de de de de de de de dede de dede L LLL LLL LL LL LLL L Yr ir Yr Yr r ir Yr oir ŵr âr ŵr ŵr ŵr yr dc dc ŵr ŵr rr dc ŵr ŵr ŵr ŵr ir ic ŵr cr dc cr oir ŵr ir Yr ŵr ŵr ŵr rr ŵr ŵr ŵr rr ŵr ŵr ir rr rr rr rr rr ŵr ŵr ŷr r ŵr ŵr 

e e de de de e rc ŵr ŵr cc dŵr جلت‎ КК ИК Ж Ж її її її її її її її Кї її К Кї її К її ŵr vir oir oir ir oir ŵr cr he ŵr ŵr Ŵr dr die oir ŵr vir oir oir ŵr ŵr câr ŵr ŵr ir ŵr dŵr ŵr ŵr rr ŵr ŵr ŵr ŵr rr dr ŵr ir ŵr ŵr ŵr ir Ŵr ŵr ir ŵr ŷr ŵr ŵr 
GENERATE TWO U(O,1) RANDOM NUMBERS 

ŵr ŵr ŷr ŵr ŷr جا‎ hee de ee de ede dede fede de dee ŵr ŷr oir ŷr ŷr oir âr e ŵr ŵr ŵr ŵr Ŵr wr ŵr rr ir ŵr r r ir dr ŵr Ŵr dc ŵr Ŵr ŵr dr ir lr ŵr dr ŵr ŵr ŵr dr ŵr ŵr le ŵr iv ŵr ŵr ŵr ŵr rc ŵr ŷr ŷr ŵr ŵr ŵr ŵr ŵr ŵr ŵr ŵr oir oir ŵr oir ŵr ŵr oir oir oir âr e âr ŵr ŵr ŵr ŵr ŵr ŵr ŵr 
RONE=GGUBFS(NSEED) 
RTWO-GGUBFS(NSEED) 
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GENERATE TWO N(0,1) RANDOM VARIABLES 
Se de de de de e de ede de dede de dede ede ee k k k k k k k k k k k k k k k k K k k k k k k k k k e ee k k k k k k k ŵr ŵr ŵr K K K R K K K k k k k K k k k k k Kk k k Kk k k k k k kk k k k k k k k le k k k K k k K k 
ZONE=SQRT(-2*ALOG(RONE))*COS(2*PI*RTWO) 
ZTWO=SQRT(-2*ALOG(RONE))*SIN(2*PI*RTWO) 
Wc ŵr de e ŵr k k k k kt fr ir fr fr fr k k k k k k k k k k k k k k Ŵr k k k k k k k k k k k k k K K k k k k k k k k k k k k kk k k Kk k k k k k k k k k k kk ŵr k k k kk ŵr kk k k k k k k k k k ŵr kk ŵr ŵr k k k k k k k K ŵr 
COMPUTE NORMAL RANDOM VARIATE WITH PARAMETERS (MEAN,SIGMA) 
FOR MEAN AND STANDARD DEVIATION 
e e de de de k k k k k k k k k k kt k k k k fk k k kk k k k k Kk k k k k k k Kk k k kk k k k ج30‎ k k k k k k k k k k k k k K k k k K k k k k k k k k k k K k k k k K k k k K k K k k k k k k kO K k K k k K k k K ŵr ŵr 
NORML=ZONE*SIGMA+MEAN 
K=l 
RETURN 
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COMPUTE NORMAL RANDOM VARIATE WITH PARAMETERS (MEAN,SIGMA**2) 
FOR MEAN AND VARIANCE 
katk dkk WWW RE WWW WR RR WR RW RW EE EERE جلا‎ dr dr dr dd od i جلا جلا لا جل جلا جلا جلا جلا‎ EEE EE EEE 
10. NORML-ZTWO*SIGMA*MEAN 
K=0 
RETURN 
END 
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IMSL ROUTINE NAME 
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COMPUTER 


LATEST REVISION 


PURPOSE 
USAGE 
ARGUMENTS — GGUBFS 
DSEED 
PRECISION /HARDWARE 


REQD. IMSL ROUTINES 


NOTATION 


COPYRIGHT 


WARRANTY 


IBM/SINGLE 
JUNE 1, 1980 


BASIC UNIFORM (0, 1)RANDOMNUMBER 
GENERATOR-FUNCTION FORM OF GGUBS 


FUNCTION GGUBFS (DSEED) 


RESULTANT DEVIATE. 

INPUT/OUTPUT DOUBLE PRECISION VARIABLE 
ASSIGNED AN INTEGER VALUE IN THE 
EXCLUSIVE RANGE (1.DO, 2147483647.DO 
DSEED IS REPLACED BY A NEW VALUE TO 
USED IN A SUBSEQUENT CALL. 


SINGLE/ALL 
NONE REQUIRED 


INFORMATION ON SPECIAL NOTATION AND 

CONVENTIONS IS AVAILABLE IN THE MANUAL 
INTRODUCTION OR THROUGH IMSL ROUTINE 
UHELP 


1978 BY IMSL, INC. ALL RIGHTS RESERVED 


IMSL WARRANTS ONLY THAT IMSL TESTING 
HAS BEEN APPLIED TO THIS CODE. NO 
OTHER WARRANTY, EXPRESSED OR IMPLIED, 
IS APPLICABLE. 


REAL FUNCTION GGUBFS (DSEED) 
SPECIFICATIONS FOR ARGUMENTS 
DOUBLE PRECISION DSEED 
SPECIFICATIONS FOR LOCAL VARIABLES 
DOUBLE PRECISION D2P31M,D2P31 
D2P31M=(2**31) - 1 
D2P31 =(2**31)(0R AN ADJUSTED VALUE ) 
DATA D2P31M/2147483647.DO/ 
DATA D2P31 /2147483648.DO/ 
FIRST EXECUTABLE STATEMENT 
DSEED = DMOD( 16807.DO DSEED,D2P31M) 
GGUBFS = DSEED / D2P31 


RETURN 
END 
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GGUNOO10 
GGUN0020 
GGUN0030 
GGUN0040 
GGUNO050 
GGUN0060 
GGUN0070 
GGUNO080 
GGUNO090 
GGUNO100 
GGUNO110 
GGUN0120 
GGUNO 1 30 
GGUN0140 
GGUNO150 
GGUNO160 
GGUNO170 
GGUNO 1 80 
GGUNO190 
GGUNO200 
GGUNO210 
GGUNO220 
GGUNO230 
GGUNO240 
GGUNO250 
GGUNO260 
GGUNO270 
GGUNO280 
GGUNO290 
GGUN0300 
GGUN0310 
GGUNO320 
GGUN0330 
GGUN0340 
GGUN0350 
GGUNO360 
GGUNO3 70 
GGUN0380 
GGUN0390 
GGUN0400 
GGUNO4 10 
GGUN0420 
GGUN0430 
GGUN0440 
GGUNO450 
GGUN0460 
GGUN04 70 
GGUNO480 
GGUN0490 
55000500 
GGUNOS 10 
GGUNO5 20 


APPENDIX B 


SIMULATION MODEL OF AIRCRAFT SIMULATOR MAINTENANCE 


LAMBDA USED FOR OFT1 FAILURE TIMES 25.80 
LAMBDA USED FOR OFT2 FAILURE TIMES 25.80 
LAMBDA USED FOR WSTÎ FAILURE TIMES 22.64 
LAMBDA USED FOR WST2 FAILURE TIMES 22.64 
LAMBDA USED FOR PTT1 FAILURE TIMES 28.05 
LAMBDA USED FOR PTT2 FAILURE TIMES 28.05 
LAMBDA USED FOR PTC1 FAILURE TIMES 
LAMBDA USED FOR PTC2 FAILURE TIMES 

MEAN SERVICE TIME FOR OFT1 REPAIR 

MEAN SERVICE TIME FOR OFT2 REPAIR 

MEAN SERVICE TIME FOR WST1 REPAIR 

MEAN SERVICE TIME FOR WST2 REPAIR 

MEAN SERVICE TIME FOR PTT1 REPAIR 

MEAN SERVICE TIME FOR PTT2 REPAIR 

MEAN SERVICE TIME FOR PTCI REPAIR 

MEAN SERVICE TIME FOR PTC2 REPAIR 
STANDARD DEVIATION FOR OFTI REPAIR 
STANDARD DEVIATION FOR OFT2 REPAIR 
STANDARD DEVIATION FOR WST1 REPAIR 
STANDARD DEVIATION FOR WST2 REPAIR 
STANDARD DEVIATION FOR PTT1 REPAIR 
STANDARD DEVIATION FOR PTT2 REPAIR 
STANDARD DEVIATION FOR PTC1 REPAIR 
STANDARD DEVIATION FOR PTC2 REPAIR 
NUMBER OF MAINTENANCE PERSONNEL 

MAXIMUM QUEUE LENGTH 

OFT1 OPERATIONAL READY RATE 

OFT2 OPERATIONAL READY RATE 

WST1 OPERATIONAL READY RATE 

WST2 OPERATIONAL READY RATE 

PTT1 OPERATIONAL READY RATE 

PTT2 OPERATIONAL READY RATE 

PTCI OPERATIONAL READY RATE 

PTC2 OPERATIONAL READY PATE 

NUMBER OF MAINTENANCE ACTIONS THAT TOOK 
MORE THAN FOUR (4) HOURS TO CORRECT 0.00 
SIMULATION RUN TIME 6000.00 
NUMBER OF MAINTENANCE ACTIONS DURING SIMULATION 550 
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SIMULATION MODEL OF AIRCRAFT SIMULATOR MAINTENANCE 


LAMBDA USED FOR OFT1 FAILURE TIMES 25.80 
LAMBDA USED FOR OFT2 FAILURE TIMES 25.80 
LAMBDA USED FOR WST1 FAILURE TIMES 22.64 
LAMBDA USED FOR WST2 FAILURE TIMES 22.64 
LAMBDA USED FOR PTT1 FAILURE TIMES 28.05 
LAMBDA USED FOR PTT2 FAILURE TIMES 28.05 
LAMBDA USED FOR PTC1 FAILURE TIMES 21.99 
LAMBDA USED FOR PTC2 FAILURE TIMES 21.99 
MEAN SERVICE TIME FOR OFT1 REPAIR 1.14 
MEAN SERVICE TIME FOR OFT2 REPAIR 2.08 
MEAN SERVICE TIME FOR WST1 REPAIR 1.46 
MEAN SERVICE TIME FOR WST2 REPAIR 1.46 
MEAN SERVICE TIME FOR PTT? REPAIR 2.23 
MEAN SERVICE TIME FOR PTT2 REPAIR 2.23 
MEAN SERVICE TIME FOR PTCI REPAIR 1.05 
MEAN SERVICE TIME FOR PTC2 REPAIR 1.05 
STANDARD DEVIATION FOR OFTI REPAIR 0.22 
STANDARD DEVIATION FOR OFT2 REPAIR 0.36 
STANDARD DEVIATION FOR WST1 REPAIR 0.43 
STANDARD DEVIATION FOR WST2 REPAIR 0.43 
STANDARD DEVIATION FOR PTT1 REPAIR 0.74 
STANDARD DEVIATION FOR PTT2 REPAIR 0.74 
STANDARD DEVIATION FOR PTC1 REPAIR 0.07 
STANDARD DEVIATION FOR PTC2 REPAIR 0.07 
NUMBER OF MAINTENANCE PERSONNEL 3 
MAXIMUM QUEUE LENGTH 1 
OFT1 OPERATIONAL READY RATE 97.% 
OFT2 OPERATIONAL READY RATE 86.5 
WST1 OPERATIONAL READY RATE 96.% 
WST2 OPERATIONAL READY RATE 95.5 
PTT1 OPERATIONAL READY RATE 94 .% 
PTT2 OPERATIONAL READY RATE 95.2% 
PTCI OPERATIONAL READY RATE 98.% 
PTC2 OPERATIONAL READY PATE 98.% 
NUMBER OF MAINTENANCE ACTIONS THAT TOOK 

MORE THAN FOUR (4) HOURS TO CORRECT 0.00 
SIMULATION RUN TIME 6000. 00 
NUMBER OF MAINTENANCE ACTIONS DURING SIMULATION 553 
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SIMULATION MODEL OF AIRCRAFT SIMULATOR MAINTENANCE 


NUMBER OF MAINTENANCE PERSONNEL 
MAXIMUM OUEUE LENGTH 


OFT1 
OFT2 
WST1 
WST2 
РТТ1 
РТТ2 


OPERAT IONAL 
OPERAT IONAL 
OPERAT IONAL 
OPERATIONAL 
OPERATIONAL 
OPERATIONAL 


READY 
READY 
READY 
READY 
READY 
READY 


RATE 
RATE 
RATE 
RATE 
RATE 
RATE 


LAMBDA USED FOR OFT1 FAILURE TIMES 25.80 
LAMBDA USED FOR OFT2 FAILURE TIMES 25.80 
LAMBDA USED FOR WST1 FAILURE TIMES 22.64 
LAMBDA USED FOR WST2 FAILURE TIMES 22.64 
LAMBDA USED FOR PTT1 FAILURE TIMES 28.05 
LAMBDA USED FOR PTT2 FAILURE TIMES 28.05 
LAMBDA USED FOR PTC1 FAILURE TIMES 21.99 
LAMBDA USED FOR PTC2 FAILURE TIMES 21.99 
MEAN SERVICE TIME FOR OFT1 REPAIR 1.14 
MEAN SERVICE TIME FOR OFT2 REPAIR 2.08 
MEAN SERVICE TIME FOR WST1 REPAIR 1.46 
MEAN SERVICE TIME FOR WST2 REPAIR 1.46 
MEAN SERVICE TIME FOR PTT1 REPAIR 2.23 
MEAN SERVICE TIME FOR PTT2 REPAIR 2.25 
MEAN SERVICE TIME FOR PTCI REPAIR 1.05 
MEAN SERVICE TIME FOR PTC2 REPAIR 1.05 
STANDARD DEVIATION FOR OFTI REPAIR 0.22 
STANDARD DEVIATION FOR OFT2 REPAIR 0.36 
STANDARD DEVIATION FOR WST1 REPAIR 0.43 
STANDARD DEVIATION FOR WST2 REPAIR 0.43 
STANDARD DEVIATION FOR PTT1 REPAIR 0.74 
STANDARD DEVIATION FOR PTT2 REPAIR 0.74 
STANDARD DEVIATION FOR PTC1 REPAIR 0.07 
STANDARD DEVIATION FOR PTC2 REPAIR 0.07 

4 

0 

7 

7 

8 

6 

4 

4 

8 


PTCI 
PTC2 


OPERATIONAL 
OPERATIONAL 


READY 
READY 


RATE 
RATE 


9 9 9 996999 
PE HL TL HL IL TL Ie بج‎ 


NUMBER OF MAINTENANCE ACTIONS THAT TOOK 
MORE THAN FOUR (4) HOURS TO CORRECT 
SIMULATION RUN TIME 6000.00 
NUMBER OF MAINTENANCE ACTIONS DURING SIMULATION 552 
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